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| | _ Halbscheibe mit Randglied. | 

Ein Spannungsgleichnis zum Problem der tragenden Linie. 
Von A. Pflüger in Hannover. | 


N Zu dem Problem der tragenden Linie der Aerodynamik wird ein Gleichnis besprochen, das durch den 
Spannungszustand einer Halbscheibe mit Randglied dargestellt wird. Der Spannungszustand besteht in einer 


glied bedingten Störung dieser konstanten Spannung. 


Be With regard to the problem of the supporting line in aerodynamics an analogy is discussed, represented 

E52” by the flow of stresses in a semi-infinite plate provided with a fish plate. The state of temsion is composed 

‘of a constant longitudinal tension parallel to the boundary line of the disk, and of the disturbance of 
the constant tension concerned, caused by the fish plate. 


_ Concernant le probleme de la ligne de soutien adrodymamique, on considere une analogie se rapportant & 
l’etat de tension d’une demi-disque avec renforcement de l’arete. L’etat de tension se compose d’une tension 
longitudinale constante, dirigee parallelement a l’ar&ie limite du demi-disque, et d’une perturbation de celle-ci, 
produite par l’element renforcateur. 

B csas3u c asponuHamnyeckoü sanayei ONOPHOÄ IHHHH BEIBONUTCH CPABHeHHe, OTHOCH- 
INeecH K HAIPAKEHHOMY COCTOAHHIO IOJLY-AHCKA C YCHIUTENBHOH HAKIANKON Ha Kkpam. Haıpa- 
3KeHHOE COCTOAHHE CAATAETCH U3 TOCTOAHHOTO IPONONBHOTO HANPAIKEHHA, TIAPAJLIEIIBHOTO 
Epalo HOJIY-AHCKA, H HAPYIIeHHSIM ITOTO HAIPAKEHHA, BEISBIBAEMOTO HAJIMYHEM HAKIAIEH. 


1. Aufgabenstellung und Bezeichnungen. 

In der Aerodynamik wird eines der wichtigsten Probleme der Tragflügeltheorie durch die 
bekannte Prandtlsche Integrodifferentialgleichung der tragenden Linie verkörpert. Auf 
formal dieselbe Gleichung führt in der Theorie der Scheiben ein besonders für den Holzbau be- 

- deutungsvolles Spannungsproblem, dessen Erläuterung Gegenstand der nachstehenden Ausfüh- 
rungen sein soll. EN DR , 

Bei vorwiegend durch Längskräfte beanspruchten Stäben, die etwa nach Bild I die Gurte 
eines hölzernen Biegeträgers sein können, wird es häufig erforderlich, laschenartige Beschläge auf- 
zuleimen, die z. B. zur Einleitung von Kräften dienen. An den Enden solcher Laschen treten bei 
Belastung des Trägers in der Regel erhebliche Spannungen auf, die sich mit der elementaren 
Balkenbiegungslehre nicht erfassen lassen, und zu deren Ermittlung besondere Untersuchungen 


Bild1. Hölzerner Biegeträger mit aufgeleimter Lasche. Bild 2. Halbscheibe mit Randglied. 


notwendig sind. Beschränkt man sich auf Verhältnisse, bei denen die Laschendicke klein ist 
gegenüber der Dicke des Stabes, auf den die Lasche geleimt ist, so kann zur Klärung des sich ein- 

stellenden Spannungszustandes folgendes Problem betrachtet werden (vgl. Bild2): 
Eine Scheibe von der Form einer Halbebene trägt ein eindimensionales Randglied von im 
allgemeinen endlicher Länge, das als völlig biegeweicher Stab aufgefaßt wird, der quer zu seiner 
Achse keinerlei Kräfte aufnehmen kann. Das im folgenden durchweg als Lasche bezeichnete Rand- 
glied erfährt also nur eine achsiale Beanspruchung, die durch Schubkräfte hervorgerufen wird, 
die zwischen Scheibe und Lasche wirken. Diese Schubkräfte hat man sich genau in der Laschen- 
achse, die zugleich die Begrenzungslinie der Halbebene ist, angreifend zu denken. Der Spannungs- 
zustand der Scheibe besteht aus einer konstanten, parallel zu ihrer Begrenzungslinie wirkenden 
Längsspannung und der Störung, die noch durch die Lasche entsteht, die die Dehnung der Rand- 
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> konstanten parallel zur Begrenzungslinie der Scheibe wirkenden Längsspannung und der durch das Rand- 


3 : schen a ni der at vom ae a 


Dicke der Scheibe und si damit übereinstimmende Breite der Lasche seie en 


- Eins gesetzt. 


Es seien die Nachktehsnter Bezeichnungen benutzt, 
‚Laschenlänge, 
IR Koordinaten a Bild 2, 


2 =.) dimensionslose Koordinaten, 
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Laschendicke für y= 0, 


Br! Elastizitätsmodul des Werkstoffes der Scheibe, 
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o Spannung o, für z=0, d.h. am Rande der Scheibe, 
E23 . Schubspannung 7 für = 0, 
o konstante Spannung in der Scheibe bei er vorhandener Lasche bzw. a 
- wert, dem sich die Spannung o, bei großer Entfernung vom Laschenbereich 
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2. Die Gleichungen für den Spannungszustand. 
In Bild 3 sind die an einem Laschenelement angreifenden, auf die Laschenbreite Eins be- 
zogenen Kräfte ee Danach ergibt sich sofort als Bedingung für das Kräftegleichgewicht 
in y-Richtung die Differentialgleichung 


d(o,t) 


+T=0. a REN 
Gt ZZ. —- ütr+a(öjt) day 
Tay . ” 5. * 


Es sei angenommen, daß die zwischen Lasche und 

- ' Scheibe befindliche Leimschicht eine vollkommene Verbin- 

Be ee Sr dung beider Teile herstellt.‘ An der Begrenzungslinie der 

Halbscheibe müssen dann die Dehnung der Faser x = 0 

der Scheibe und die Dehnung der Lasche miteinander übereinstimmen. Daraus folgt als weitere 
Beziehung 


2.05 N 
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Der in dei Halbscheibe auftretende Spannungszustand wird am nee in folgender 
Weise quellenmäßig dargestellt. Eine im Punkte x = 0, y= y’in Richtung der negativen y-Achse 
auf die Scheibe wirkende Einzelkraft P erzeugt bekanntlich die Spannungen!) 


er ey—y) 5 BES Blna9) are ey ie 
ey) 7°. me ” [8 +(y—y')2]? 
Die durch die Randschubkräfte 7 in der Scheibe hervorgerufenen Spannungen ergeben sich nun, 
wenn man bei Benutzung von y’ als Zwischenvariabler P=7dy’ setzt und über den ganzen 


Laschenbereich integriert. Berücksichtigt man noch die konstante Spannung o, so erhält man 
schließlich den gesamten Spannungszustand der Scheibe zu 


a 


!) Vgl. z.B. Föppl: Drang und Zwang. Bd.1. München und Berlin 1941, S.295, — Eine ausführliche Dar-, 


ar der Integralgleichungen des ebenen Spannungszustandes findet sich bei Weinel: Z. angew. Math. Mech. a 
(1931), S. 349. 


Vatt (y- =y ya? 


7 
Danach Tassen sich die Spannungen für jeden Punkt der Scheibe Denen wenn T = bekannt i ist. 


1 v laßt sich nach (1) aus 07, ermitteln; co; wiederum ergibt sich ne 


U Aus = folgt für = = für die EL co 


# a Talchenbereich. d.h. für. —— 2 Mi +2 ergibt Sich dabei ein uneigentliches Integral wie 


man sich bei Ausführung des eg x2—-0 überzeugen. kann, erhält man durch den - 


Cauchyschen Bew des Integrals den ENDEN Spannungszustand. Es ist also zu schreiben 
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Wird But Met genauere Schreibweise im n folgenden wieder verzichtet, so ergibt sich unter Benktzung 
von (1) und (2) aus (4) 
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Man erhält also eine dla für die Veränderliche ozt. Diese Gleichung 
stimmt nun bemerkenswerterweise formal mit der in der Aerodynamik des Tragflügels vorkommen- 
den Prandtlschen Gleichung a ren Linie überein. Schreibt man nämlich ur letztere 


in der Form 2) 


1 
Gt= — Ogoo ı 
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z *) vgl. z.B. Prandtl-Tietjens: Hydro- und Aerodynamik Bd.2. Berlin 1931, S. 211. 
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Flügeltiefe entsprechend der 


b - Flügelspannweite entsprechend der Laschenlänge b, ES N EN EA 
Ca örtlicher Auftriebsbeiwert entsprechend der Laschenspannung u 0... 
&  Anstellwinkel entsprechend 1% ARE = N 2 x i ä : & u 
cs _ Profilkonstante (Ableitung von c, nach a für unendlich große Flügelstreckung) 
entsprechend 16 Er es IF < RER Ei 


‘Die Gestalt der Lasche entspricht also direkt dem Flügelumriß, die Zirkulationsverteilung 
der Laschenkraft und die induzierte Geschwindigkeit der durch die Lasche in der Scheibe hervor- 
gerufenen Spannungsstörung. Die Analogie zwischen den beiden Problemen dürfte damit beson- 
ders anschaulich sein. Sie ist jedoch nicht nur deswegen von Interesse, weil damit ein Spannungs- 
gleichnis für das in der Tragflügeltheorie so wichtige Problem der tragenden Linie gefunden ist, 


a 


sondern vor allem auch im Hinblick auf die eingangs gestellte Spannungsaufgabe selbst. Denn die 


nicht unerheblichen Schwierigkeiten, die die Lösung der Gleichung (5) hauptsächlich infolge der 
auftretenden Singularität bereitet, erledigen sich jetzt dadurch, daß die Erkenntnisse und Methoden, 
die in der Aerodynamik bei Behandlung der Prandtlschen Gleichung gewonnen sind, hier ohne 


weiteres verwendet werden können. Mit Rücksicht auf die Spannungsaufgabe wird hier nur ein 


konstantes Verhältnis & und eine konstante Spannung o betrachtet. Das durch (5) zum Aus- 
druck kommende Gleichnis gilt natürlich auch für einen Flügel mit veränderlichem c;. und ver- 
änderlichem Anstellwinkel. Man braucht dazu nur E, als von y abhängig anzunehmen und sich 
ferner die Halbscheibe geeignet so belastet zu denken, daß für die Scheibe ohne Lasche die Rand- 
spannung dem gewünschten Verlauf von «& entspricht. °®) 
Für die weitere Untersuchung wird zweckmäßig (5) noch etwas umgeformt: 
+1 

| bu ee Bor 

OR nb E)Jdy\oh, a 
—1 


oder mit den in Abschnitt 1 angegebenen abkürzenden Bezeichnungen 


| ra 
i-rlı—al i | RE a 
ae a 
In manchen Fällen, insbesondere dann, wenn es sich in der Hauptsache nur um die Klärung 
des Spannungsverlaufes an den Laschenenden handelt, ist es wünschenswert, eine ‚‚halbunend- 
lich lange‘“ Lasche zu betrachten. Hierbei wird es notwendig, statt 7 die neue Koordinate 


a 

2 a | 
= Er RL d.h. den durch i, dividierten Abstand vom Laschenende, und statt A den neuen Para- 

0 

= 2E 7 b . | ; 

meter A = TE art zu benutzen. Nach Ausführung des Grenzüberganges b5—> oo ergibt sich _ 
[1] x “ 
dann die Gleichung 
din dt: | 
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3. Die elliptische Lasche als Beispiel. 

Über die Eigentümlichkeiten des Spannungsverlaufes in der Lasche und der Scheibe erhält 
man am besten einen Überblick, wenn man ein möglichst einfaches Beispiel für eine Lösung 
von (6) betrachtet. In der Tragflügeltheorie wird die bekannteste und auch praktisch wichtigste 
Lösung der Prandtlschen Gleichung durch eine elliptische Auftriebsverteilung mit elliptischem 
Flügelumriß und konstanter induzierter Geschwindigkeit dargestellt. Es liegt daher nahe, auch 
hier zunächst einen elliptischen Verlauf der Laschendicke zu betrachten. Wird der Ansatz 


D= LY1—n? 


°’) Im Hinblick auf das Gleichnis sei noch erwähnt, daß ein von F. Ziller I 
nbli ‚hnt, R ‚„ Ing.-Arch. 11 (1940), S. 
angegebenes Variationsproblem zur Prandtlschen Gleichung genau mit dem Casti gliano a ne a 
Minimum der ‚Formänderungsarbeit für das Spannungsproblem übereinstimmt, wo durch eine anschauliche Deutun 
der von F. Ziller rein formal Aufgestellten Beziehung möglich wird. F 
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gemacht, wobei Z, der noch zu bestimmende Wert von L für n=0ist, so wird fürn? <1 er - 


u? Aa Anz a a nähe > a a Rt 
Ri + N e r u. 
er 0 


y F n—E. rl i 
a { dn' { n' 
% Er . —1ı Yln a 1 +8 Te BE 55 
Für „=0ist L=L,und T=1; daraus folgt L,= Erle Im’ Die praktisch interessierenden Größen Be 


sind die Laschendicke und die Spannungen o und 7; sie werden am besten in der dimensionslosen 
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Form es ud angegeben. Man erhält 
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| ee 
Ta a 
DEE a I 
| Fr = al.) 1 i 
: | U u a A 
®..,0 dy 2 dn 22 1+Ar Yı m } 
I Pay ee —n?. NE WR We 3 
; % lerrT, Be | a 
| Außerhalb des Laschenbereiches für 7? >1 muß = Osein. Für die Randspannung o ergibt sich 

ee +1 | 
i Er 7 .d L Pe F 

o ee Bet | Be 


Das Ergebnis der Rechnung ist für A= 1 in Bild 4a — c aufgetragen. Im Laschenbereich 3 

ist den bekannten Ergebnissen der Tragflügeltheorie entsprechend der Verlauf der Laschendicke 5 
elliptisch und die Randspannung 5 konstant. An den Laschenenden wird o sprunghaft unendlich a 
= groß und nähert sich dann mit weiterer Entfernung vom Laschenende dem konstanten Spannungs- 


AN 


a 


Bild4a. Längsspannungen o,, für eine elliptische Lasche. 


zustand o, der für die ungestörte Scheibe gilt. Das Verhalten von o, in &-Richtung ist dadurch 
‘gekennzeichnet, daß sich ebenfalls die Randspannung mit größer werdendem & überall dem Werte o 
nähert. Die Längsspannung o, verläuft für alle Schnitte 7) = konst. derart, daß sie am Scheiben- 
rand voraussetzungsgemäß Null ist, mit wachsendem & ein Maximum erreicht und schließlich 
wieder auf Null abklingt. Die Größe von o, ist für die ganze Scheibe praktisch bedeutungslos. Die 


dLl 
Schubspannung r ist am Scheibenrand der Ableitung 3 proportional und wird. dementsprechend 


an den Laschenenden unendlich groß. Im übrigen nimmt 7 vom Randwert 7 ausgehend mit der 
Entfernung vom Laschenbereich auf Null ab. Bemerkenswert ist, daß r für Schnitte n = konst. 


DIR See, 


SER ERTRN Pflüger, Halbscheibe mit Randgli 
- das Vorzeichen wechselt und ein Minimum bzw. Maximum aufweist. Bei dem vorliegenden Bei 


en “++ £_1.0 die Schubspannung für alle n schon < 0,05 0, also praktischim 
u Se Der Spaßnunpszuktänld von Bild4giltdmnach 
Br. 


7 = 


Br 3 . Y ri Aacı 

Br Vergleich zur Längsspannung o vernachlässig 5 Aa Re. 3 ; : Re 
mit hinreichender Genauigkeit auch für eine unendlich lange streifenförmige Scheibe, die au Bu 
2. beiden Rändern symmetrisch eine Lasche trägt, wenn die Streifenbreite > 5 ist. er. 
2 = ; ia £ > £ $ 
z h ci * i Ren 
B. I i 
ra 

Es 

Be. » 
£ Er. 
A 


Bild4c. Schubspannungen , für eine elliptische Lasche. 


Die wesentliche Erkenntnis, die aus den Bildern 4folgt, ist die, daß die aufgeleimte Lascheden 
Spannungsverlauf in der Scheibe sehr ungünstig beeinflußt. Im Laschenbereich wird zwar die 
Längsspannung o vermindert; dafür treten jedoch Schubspannungen auf, die am Laschenende un- 
endlich groß werden. Ebenfalls geht die außerhalb des Laschenbereiches gültige — in der Aero- 
dynamik bedeutungslose — Lösung von o für 7?—= 1nach unendlich. Die Laschenenden sind also 
Stellen besonders ungünstigen Spannungsverlaufes und verhalten sich in ihrer Wirkung etwa wie 
Spitzkerben. Es dürfte daher von besonderem Interesse sein, den Spannungszustand am Laschen- 
ende näher zu untersuchen und festzustellen, unter welchen Umständen unendlich große Span- 
nungen vermieden werden. können. 


4. Der Spannungsverlauf am Laschenende. 
Es sei zunächst die Schubspannung 7 betrachtet. Für diese gilt die Beziehung 
- nm ,dL mA (&s do 
T=— = — [07 t 
er A. 


aus der sich sofort entnehmen läßt, daß bei dem Beispiel der elliptischen Lasche 7 deswegen un- 
endlich groß wird, weil am Laschenende die Tangente des Laschenumrisses auf der n-Achse senk- 
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recht steht und infolgedessen —— — » ist. Diese Erkenntnis dürfte ohne weiteres einleuchtend 


d 


sein, daman es gewohnt ist, bei zu plötzlichen Querschnittsänderungen von Baugliedern theo- 


retisch unendlich große Spannungen zu erhalten. Umgekehrt ist zu vermuten, daß man end- 


zum Ende hin allmählich auslaufen läßt. 


liche Schubspannungen erhält und darüber hinaus 0 endlich bleibt, wenn man die Laschendicke 


Um über diese Vermutung Klarheit zu schaffen, werden am zweckmäßigsten die Verhältnisse Ba 


an einer halbunendlich langen Lasche betrachtet, da es ja nur auf den Spannungszustand am 


Laschenende ankommt. Bei Behandlung von Gl. (7) ist es unbequem, daß sich das Integrations- 
"Intervall von Null bis Unendlich erstreckt. Eine Reduktion auf das Intervall von Null bis Eins 


wird durch Einführung einer neuen Veränderlichen 


A | | 
WERE Bar? ’ 


möglich, wobei a einen noch passend zu wählenden Parameter darstellt. Man erhält nach Durch- 


führung der Substitution aus (7) . 


I 5 
L=T|1-aı u) = dw] 
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. Eine praktisch wesentliche Erkenntnis läßt sich nun schon gewinnen, wenn man das einfache 
Beispiel : 


betrachtet. Es ergibt sich dann für£>0 aus (8) 


zulmail — u) f ne | 


1—u 
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Be 
Ts Z anı R 
oe "de (tat 
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‘und außerhalb des Laschenbereiches für & <O 


Bild5. Dicken- und Spannungsverlauf für eine halb- 
In(- 4 2] unendlich lange Lasche mit endlicher Schubspannung für 
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Der Verlauf der Laschendicke und der Spannungen in Abhängigkeit von £ istfüri=1 
und «=1 in Bild 5 aufgetragen. Die Schubspannung bleibt jetzt überall endlich, Die Rand- 
spannung o wird jedoch am Laschenende wiederum unendlich groß, und zwar geht sowohl die 
innerhalb, als auch die außerhalb des Laschenbereiches gültige Lösung über alle Grenzen. Be- 
sonders bemerkenswert ist, daß am Laschenende wegen o —= w nach der Beziehung 


+ 


| dL_ do 
RT NE Te 
I=10=, N 
SE — (Q wird. Man erhält somit das überraschende Ergebnis, daß sich die unendlich großen Span- 


nungen trotz weitgehender Ausrundung des Überganges zwischen Lasche und Scheibe ergeben. 
Es ist also nicht möglich, ein langsames Auslaufen der Laschendicke als ausreichende Vorschrift 
für eine zweckmäßige konstruktive Ausbildung des Laschenendes anzusehen. Der Grund für das 
Verhalten von oist in der Unstetigkeit der Randschubspannung 7 zu suchen, die am Laschenende 


von dem Wert 4 plötzlich auf den außerhalb der Lasche gültigen Wert Null abfällt. Man kann 
nämlich endliche Spannungen erhalten, wenn man die Funktion Z so wählt, daß auch ihre erste 
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% plane und damit die Randschubsp 
Ben Runge ergibt sich auf ee a 


Es sei vorausgesetzt, daß sich. hin im Bereich 0 ” us 
von der Form ARTE VERA N ‚= s ER Be 


Bi, OR läßt. Der ns betrachtete Fall, daß T un damit 4 Tas In x 
endlich groß werden, wird also im folgenden Be ss ” folgt dann ee 
ee A w 3 a - Enahe Fk 
GE ai wit an Kosw IE un = ‚dw. er ie 
Bei der Tnteeratian wird aus dem n-ten Glied der im Integranden stehenden Reihe R 
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nahme des letzten Gliedes Ausdrücke liefern, die am Laschenende bei u = 0 für alle in Frage 
kommenden n endliche Werte annehmen. Das letzte Glied dagegen bedingt ein Integral 


1 2 um Er hr Ca “. h ; x en oh 
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.  - Man sieht nun sofort, daß sämtliche Glieder in der Kamer unter dem Integralzeichen mit A © ER 
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das für u = 0 verschwindet, wenn n list, jedoch unendlich groß wird, wenn n = ist. Sollen 


dL 
also unendlich große Randspannungen 0 5 vermieden werden, so muß in der Reihe für Zu das Glied. 
. mit c, verschwinden oder, da ” 


Bi Sen: U 
Er, Per n T = —— — — re 
E: T | on) TE onAa(l— u) Fu 
Br ist, mit anderen Worten: Am Laschenende muß die Schubspannung 7 = 0 sein. Unter den an- 


of gegebenen Voraussetzungen stellt damit das Verschwinden der Randschubspannung am Laschen- _ 
Be... ende eine notwendige und hinreichende Bedingung für einen endlichen Spannungszustand dar. 
we Aus dem Vorhergehenden ergibt sich nun direkt, daß bei der halbunendlich langen Lasche 
=. z.B. durch den Ansatz 
2 L=0,W +0, W...+Onur | 


e-. | unendlich große Spannungsspitzen vermieden werden. Auch bei endlicher Laschenlänge lassen 
sich entsprechende Ansätze leicht finden, wenn man die Forderung T= 0 für das Laschenende 
erfüllt. Z.B. führt bei Laschen, deren Dickenverlauf zur x-Achse symmetrisch ist, der Ansatz 


= L=[0, +0, P+ Om... +0, —( 20) m], 


y 
® wobei für m beliebig die Zahlen 2, 4, 6... gesetzt werden können, zum Ziel. ‘Als Beispiel für - 
m Laschenformen mit endlichen Spannungen. sind für die halbunendlich lange Lasche in Bild 6 eine 
. Reihe von Kurven zusammengestellt, die sich auf Grund des einfachen Ansatzes Z = u? ergeben 
und für a4 = 1,6 gelten. Bemerkenswert ist dabei, daß auch hier die Tangente des Laschen- 
umrisses für © = 0 mit der ö-Achse zusammenfällt. Hierzu läßt sich auf Grund der Beziehung 
er -dL ad (gr do ) 


ron TE page! 


sagen, daß bei Ba en Sonderfälle, daß am Laschenende # = 0 ae n Te) ist, die Forde- 


\ rungT=(0 nur Er 7 m = 0 erfüllt werden kann, d.h. zur Vermeidung unendlich großer Span- 


nungen ist es im allgemeinen notwendig, die Laschendicke mit der Ableitung Null auslaufen zu 


lassen. Daß diese Bedingung jedoch keineswegs hinreichend ist; war bereits oben an dem Beispiel 
L= u gezeigt worden. 
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er - Wenn es auch keine Schwierigkeiten bereitet, Laschenformen mit endlicher Spannungs 
verteilung zu finden, so wird man doch in der Praxis häufig auf die Anwendung einer derartigen 
A ‚Form verzichten, da die.wirkliche Ausführung von Laschen, wie sie etwa durch Bild 6 gegeben sind, 


und insbesondere die Einhaltung der Bedingung = 0 am Laschenende sich aus Fertigungs- RR: 


252 220 75 10 5 -5 


. Bild 6. Dicken- und Spannungsverlauf für halbunendlich 
lange Laschen mit endlichen Spannungen für Z =u® und 
5 aiA=1,6. 


gründen vielfach verbieten dürfte. In diesen Fällen muß man dann unendliche Spannungen in 
Kauf nehmen. In Wirklichkeit gilt dabei natürlich die Theorie nur bis zu einer gewissen Umgebung 
der singulären Stelle. Die meist geringe Größe dieses Bereiches und die in ihm auftretenden Span- 
nungen müssen durch Versuche geklärt werden, wie es von ähnlichen Problemen der Kerbspan- 
nungslehre her bekannt ist. 


5. Zusammenfassung. 


Wird eine Halbscheibe, die parallel zu ihrer Begrenzungslinie durch eine überall konstante 
Längsspannung beansprucht wird, mit einem biegeweichen Randglied (Lasche) versehen, so 
entsteht dadurch eine mit der Entfernung vom Laschenbereich abklingende Störung des kon- 
stanten Spannungszustandes. Die Berechnung dieser Störung führt auf eine Integrodifferential- 


‚gleichung für die Laschenkraft, die formal mit der Prandtlschen Gleichung der tragenden 
Linie übereinstimmt. Es wird daher ohne weiteres die Übertragung der Methoden und Ergebnisse _ 


der Aerodynamik auf das Spannungsproblem möglich. Lediglich der Spannungszustand am 
Laschenende erfordert einige besondere Betrachtungen. Es zeigt sich, daß hier die Spannungen 
in der Regel theoretisch unendlich groß werden, daß also praktisch eine Art Kerbwirkung auf- 
tritt. Es läßt sich jedoch auch eine bestimmte Klasse von Laschenformen mit endlichen Span- 
nungen finden, die in der Hauptsache dadurch gekennzeichnet sind, daß am Laschenende die 
zwischen Scheibe und Lasche übertragene Schubspannung verschwindet. Ein langsames Aus- 
laufen der Laschendicke so, daß ihre Ableitung am Laschenende zu Null wird, ist dabei keine 
hinreichende, sondern nur eine im allgemeinen, d.h. bei Ausschluß gewisser Sonderfälle, not- 


wendige Bedingung für die Vermeidung unendlich großer Spannungen. 


Abgeschlossen am 1.4. 1944. 


Über die Unstg der zweiten Abl rc 
 Schwerepotentials. an Diskontinuitätsflächen der | 
Von Hans Driel in Berlin. | & 


Es wird eine Formel entwickelt, welche die an einer Diskontinwitätsfläche der Dichte. ee =, 
malentrichtung (n) auftretenden Sprungwerte der zweiten Ableitungen des Schwerepotentials nach zwei b 
liebigen Richtungen (r,s) als Funktion des Dichtesprungs und der Richtungen (r,s; n) darstellt. Er 

In the case that discontinuous density is arising at a superficies with Ihe normal.direction (n), a formula 
is developped giving the leaping values of the second derivatives of the potential of gravity in two o,anbilrarg irect- 
ions (r, 3) as functions of the leap of density and the directions (r, s;n). 3% 

On deduit une formule, qui, suivant deux directions arbilraires (r, 8) et en fonction de variations brusques Ru 

- de la densite, ainsi que des directions (r, s; n), donne les valeurs correspondentes de variations discontinues 
de la deuxieme derivde du potentiel du champ de BLORHLNIER se prodwisant le ee de la surface de a 

. de la densite avec la direction normale (n). 

Beisonurca dopMyna, KOTopaA, No MByM IPONSBONBHEIM Hanpaprennam (r, 5), zaer 
SHAYEHHA Pe3KUX CKAYKOB BTOPOH IIPOH3BONHOÄ TIOTEHIHAIA CHIIEI TABKECTH, HAaÖMONaeMEIX S 

‘ BAONb  IOBEPXHOCTH HIPePBIBUCTOCTH IINIOTHOCTH © HOPMANIDHEIM BR (n), zack = 5 
PYHEIMIO PesKUX CKAUKÖB INIOTHOCTH U Hanpapıennä er 5; n).' \ - | 


ns A 2 ey 5 =, + 
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Füie Gleichung einer Fläche, a an der sich die Dichte o(z, y,z) der Materie unstetig ändere (Dis- 
kontinuitätsfläche nullter Ordnung), bezogen auf ein mit der rotierenden Erde fix verbundenes, 
orthogonales kartesisches Koordinatensystem (x, y, 2), dessen z-Achse mit der Rotationsachse 
‘der Erde zusammenfällt. Die beiden Seiten von F mögen durch die Indizes 1 und 2 unter- 
schieden und die Sprungwerte.einer physikalischen Größe y(x, y,2) beim DunDeE durch # 

‘ mittels geschweifter Klammern {4} en werden: 


ee ee 
Die Potentialtheorie zeigt, daß. 2 Schwereotentil Rt 3 
D (x, y,2 [| e a4: ee RE ee 2), 


Be — Gravitationskonstante, = Winkelgeschwindigkeit der Erdrotation um die 2-Achse), 
welches in den Räumen außerhalb F der Poissonschen Gleichung 


5 AB= ter £ +) 9= ARRO—-AON . E RE dN 
genügt, mit seinen ersten Ableitungen, nr Kraftkomponenten 
x=- 2-9, 7=-—.--9, en, En | 
beim PETE des Aufpunktes durch F stetig bleibt: 
ID SO ee &; (6); 
ee = 25 


während die zweiten Ableitungen (B,,, D,, —= ©®,, usw.) Unstetigkeiten aufweisen. it. rs] ) 


Zwecks Berechnung dieser Unstetigkeiten, d.h. der Sprungwerte {®,,}, 2 Pu} 
usw. betrachten wir zwei infinitesimal benachbarte Be P(z,y,2) und P(e+de, y+dy, 


2 -+dz) auf der Fläche F(x, y,z) = 0. Der Vektor En = (dx, dy, dz) liegt in der die Fläche F 
in P berührenden Tangentialebene, d. . es gilt die Gleichung: 


or or 
ir + +, u + Zde=M, de +F,dy+R,d=0.......(8) 


“ ” ” . “ * . er . 
im übrigen ist die Orientierung von ds in dieser Tangentialebene beliebig. 


Für eine an der Fläche F stetige Größe, z.B. X = —©,, muß gelten: 
> x 


; . (dX),= (iR). 
unter dX die Änderung von X beim Übergang von P(&, y,2) zu P(e-+de, y+ dy, 2+.dz) ver- 


standen, also 
er), l)a un. 


an HH), e= 


Er 


G1. 9 daß der, are 2 Gradienten v von “ : 


age, ie” 
ER RE = a 


chenormak n von F im Punkte Pia Y, 2) fällt, deren. Richtung durch ER 2 


ie 


Ne SR T Fa) ae: in are. a 
RER FE ‚cos m: a cos (m, Y); 2 (n, = = Fr? 55 7 2% NET = Re 
Br x gegeben 2 Es muß a das Vektorprodukt [grad r, ed 2} nen 


ter ER Ps eradR. {grad X}] = De A BON ER 
Zu weiche Bedingung en Gleichungssystem äquivalent ist: en | RN 
| ae BR JAXYT: 
| . .taR OR - “ aa 
3 . Fi aa A Fa as BR ER ie 
A a 9 2. =) = 
das aber nur zwei unabhängige Gleichungen darstellt. Ber 
; . Analoge ET, angewandt auf no Kraftkomponenten Y zung Z, führen auf Mr Be: 


- 


Be Gleichungen: | Ba 
Bi ea = | ae a a 
B:: 2 Ä Bee NET F, - \0e Nr oy) er 
ee 2 ee an, = KON 

2 KarıT y del 02 

EB ni 

Bloom ade) Br (5 
——— = r Er (17) 

& F, a) = F, a 

Kir: ; °y 02 02 
er Er Aus den Systemen (15), (16), (17) folgen dann nachstehende Beziehungen zwischen 
den ersten Ableitungen der Kraftkomponenten an der Diskontinuitätsfläche: e 
Es ER, efor) [er or far) » EN 
E. A as. 

Be ee a. 

1 0% 02) \öxJ - 
Eur; | ir a = 2 nn BE SE 20% 
E 5 0Y E79] 02 au) j 
3 - die nach Einführung des Schwerepotentials ® in die Form | 
3 SET EL BASE (21), 
# DD ED NAD ELBE CR (22), 20 = 
{ RE ' {D,,} {® 2} = {9,2}? = = {9}? ee een Se RB te (23) Er 
E Abefechen (Weitere drei Beziehungen am Schluß der Arbeit.) : 
Er Nehmen wir die aus der Poissonschen Gl. (4) folgende Beziehung 
E : Bar +@gstnal. : een. .. (24) 
En: hinzu, : so resultieren durch paarweise Addition der Gl. (21), (22), Se die Ausdrücke: 
j Anl Dr =D HOP FR en » (25), 
4 in a 
E:: aD Br FI HD Se (27): 


welche zunächst zeigen, daß die Sprungwerte {B,.}, {9,, „‚ und {®,} das gleiche Vorzeichen, 
nämlich das Vorzeichen des Dichtesprungs aufweisen müssen. 
Weiter ergibt sich aus den Gl. (13) und (15), (16), (17) mit Rücksicht auf Gl. (25), (26), (27): 


a RE er er . (28), 
ne I LS a (29), 
nn IE SEHE re ee (30), 


ER wa 


oA, die Dre G% Be (30) durch Kom ination der Gl. ( 
"Bildet man aus Gl. (13) die Produkte cos (n, ©) 

BE (n, 9) - cos (n, 2), so erhält man durch eine an R | 

een zweiten Ableitungen des Potentials die Beziehungen 


- (= Az „fo: cos (n, 2%) ) cos En. | 


DO = :{8,}= An {o}cos (n, x) cos mi) 
{9,.} = {9,} = 47% {o} cos (n, y) cos (n,2) - 


& welche nn alich aus den Systemen (21), (22), (23) und (28), (0) ve res ie 


Mit E Bezeichnungen = %, BR FD 2 a: sich AR die, Ei (28 


e and a a (33) in. BE | a Er 


Kr elf 328 57 
| bi : a” | = arten. I im, 2) 


;,k=1,2,3) zusammenziehen. | a 


Es seien nun r und s zwei beliebige, nicht notwendig orthogonale Richtungen. Wegen B* 
a 


j 6) 
Ir 2 cos, er; E FR 


N 3. 3 3 (s a Ge > 
Te ee > | k 02% A 
| 2 
ne =® 5 S Elan (s, ©) nn (35) Re 
| Ara 08or a (no ano = 
und somit Sa a RE REN 
BE) RX) 3 ’ 3 Be: H ER 2 z 
PP TEE IE) Esser 
Die Substitution von Gl. (34) in a führt dann 2 die Gleichung en j 
Ep Scost (n, &;) cos (r, ©;) * cos (N, %,) cos (& ee (BD = 
die sich wegen | } 
3 x 
cos (n, r) = Si cos (N, &;) cos (r, &;) is ErEE 
und 
3 3 
cos (n, s) = PN cos (N, %;) cos (8, %;) 
auf 


7 Ö DER. : BE 
a = al {0} cos (n, r) cos (n, s) ee 


reduziert, welche Gleichung die Sprungwerte der zweiten Ableitungen des Schwere- 
potentials nach zwei beliebigen Richtungen (r, s) an einer Diskontinuitätsfläche 


der Dichte mit der Normalenrichtung (n) als Funktion des Dichtesprungs {g} 
und der Richtungen (fr, s; n) darstellt. 


Da das Schwerepotential ® sich additiv aus dem Gravitationspotential (Attraktions- 
potential) ® und dem Zentrifugalkraftpotential ® zusammensetzt und die zweiten Ableitungen 


7 
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Ar 


a a da cha date: ZU aan 4 Haan Lu) Län N Das al En u 2 m D 


iz Be Per Di 
Ü Br N ER 2) AP u 3a KH 

ei aa pi Aue 
a War L 1% 


x RE ae a Kr Be et f 
Unstetigkeiten der zweiten Ableitungen des Schwerepotentia 


le ie a EN 
u. RT En Fe 
ED. 7 EISEN N A Be 


Über die er | 


or ds) 


von ® durchwegs stetig sind: en = 0, gilt die Beziehung Gl. (38) auch für das Gravitations- 


> potential ® allein. 2 


Die vorstehend abgeleiteten Gleichungen ermöglichen ferner die Aufstellung einer ein- 


fachen Formel für den Sprung der-mittleren Krümmung 
Be ; VER fe 1 1 
nn K ee 
3 n | : R, RR. 
einer Äquipotentialfläche ® = const. beim Schnitt mit einer Diskontinuitätsfläche der Dichte. 
Es sei P ein Punkt der Schnittkurve der Flächen ©(z, y, 2) = const. und F(x, y,2) = 0, N die 
Normale von ® in P (positiv in Richtung zunehmender ®-Werte gezählt), die wir zugleich 
zur +2-Achse eines orthogonalen kartesischen Koordinatensystem machen, dessen x, y-Ebene 
also Tangentialebene der Fläche ® = const. im Punkte P ist. Die Normale von F(x, y,2) = 0 


: ur h 0 ne 

| in P werde wie bisher mit n bezeichnet. Esist nun, unter g = 22 == = die Schwerebeschleuni- Be; 
_ gung verstanden: DE 
1/28 - 
Eat 
also, da g beim Durchgang durch F stetig bleibt: 
1 (28 ®® L 9107 z 
IK} = Rn a 4) —= re Cru ge 3 Ka ale (77) (39). 


Nun ist weiter: 


_ 608? (n, x) + cos? (n, y) = 1 — cos? (n, 2) = sin? (n, N), 
womit man aus Gl. (39) die einfache Gleichung 


Kg N FR AN 


zur Berechnung des Sprungwertes der mittleren Krümmung einer Äquipotentialfläche (Geoid- 
‚ fläche) beim Schnitt derselben mit einer Diskontinuitätsfläche der Dichte erhält. 
Die vorstehend abgeleiteten Formeln sind u.a. für Fragen der physikalischen Geodäsie 
und Geophysik von Bedeutung, worauf aber an dieser Stelle nicht eingegangen werden soll. 


Aber zur Erleichterung der Übersicht über die vorstehenden Rechnungen sei noch Folgendes 


bemerkt: Die Berechnung der Sprungwerte der zweiten Ableitungen {®,}—= {®,,} des Schwere- 
potentials nach zwei beliebigen Richtungen (r,s) an einer Diskontinuitätsfläche der Dichte 
beruht auf der Bestimmung der sechs Sprungwerte {®,}= {B,} (Ü, k= x, y, 2), wozu sechs 
unabhängige Gleichungen benötigt werden, die folgendermaßen gewonnen wurden: 

Jedes der drei als Folge der Stetigkeit. der Kraftkomponenten (X, Y,Z) an der Diskonti- 


nuitätsfläche resultierenden Gleichungssysteme (15), (16), (17) enthält zwei unabhängige. 


Gleichungen, aus denen zusammen sechs Beziehungen zwischen den sechs Sprungwerten der 
zweiten Ableitungen {®,,}={®,;} folgen. Von drei dieser Beziehungen haben wir in dem 
System (18), (19), (20) bzw. (21), (22), (23) Gebrauch gemacht, indem diese drei Beziehungen 
zusammen mit dem System (13), welches zwei unabhängige Gleichungen enthält und mittels 
.der Gl. (24) (zusammen also sechs unabhängige Gleichungen) gestatten, die-sechs Sprungwerte 
{8,,} = {9,,} durch die Normalenrichtung (r) der Diskontinuitätsfläche und durch den Dichte- 
sprung {e} an derselben gemäß Gl. (34) auszudrücken. Die restlichen drei Beziehungen zwischen 
den {B;,} (uk = x, y, 2) lauten: 73 
{D..} = {Da} {Dal {D,.} > 
Dun} = 1Dya  1Dı I Bat » 
19.) = Bar Da I {Bau » 
und werden durch die Lösungen Gl. (28), (29), (30) und (31), (32), (33) erfüllt. 
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n. Überschallströmungen pflanzen sich Unstetigkeiten längs der Charakteristiken fort. Die Art der Fort- R 
en wird für an an des ebenen und rotationssymmetrischen LT Überschall- E 
und Schallgeschwindigkeit untersucht. Beim Strahl mit Schallgeschwindigkeit gibt n eine a: 
deren einer die Entwicklung einer Überschallströmung aus einem Parallelstrahl mit Schall w ren @ gr 
Entstehung eines Parallelstrahls mit Schall aus einer ae Zn ae of Pe 3 ’ 

In suprasonic flows, discontinuities are transmitted along the characterıslıcs. Lhe kin EN I 
(R einien in ik special case of the plane and, the. rotary-symmetric parallel jet moving with suprasonic 
and sound velocities. Im the last case, there are two types of solutions the one of which represents a ee 
current developping out of a parallel jet moving with sound velocity, while.the other represents a parallel je 2 
moving with sound velocity developping from a subsonic flow. 4 : + 
i Dans les propagations & ultra-son les discontinuites se propagent le long des a La Se 
u. propagation est analisee poar le cas particulier de rayons paralleles ‚plans ou rotatif-sym triques, er : : ar, S: 
-  detransmission d’ulira-son et de son. Le rayon ü vitesse de transmission de son condwit a 2 types de soluti “ = 
er dont l’une deduit ume onde d’ultra-son d’um rayon parall&le de son, tandisque l’autre montre la formation 
celui-ci en partant d’une onde de sous-son. : Er 

B TeyeHuusx YAIBTpa-sByKa - TePEePEIBEI PACIPocTpaHamTca BAONB XAPaKTepHCTHK. Porn . 
 PacupocTpaHeHunad PACcMaTpuBaeTca MIA YACTHOTO CAIYYAA INIOCKOTO MHIIH POTAIHOHHO- 
CHMMETPHYKOTO HAPAIIENBHOTO Iy4a CO CKOPOCTEO YJABTPa-3ByKka U 3ByKa. Ina ayya co 
CKOPOCTBO 3Byka HMeIMTcA Ba Tuna pemeHnä: ONUH US HHX IIOKaskIBaeT O0PaBOBaHHe = 
TeyeHHA YIETPA-3BYKA U3 TAPAILTEIBHOTO Iy4a CO CKOPOCTBO 3BYKA, a BTOPOf-o0pasoBaHHe 
HocAeNHero H3 TEYeHHA CO CKOPOCTBP HHIKe CKOPOCTH 3BYKA. 


1. Einieitung. 

Bei der Konstruktion von Überschalldüsen pflegt man die Strömung im engsten Querschnitt 
als Parallelströmung mit Schallgeschwindigkeit aufzufassen. Bei der nachfolgenden Erweiterung 
steigt die Geschwindigkeit und es entwickelt sich ein Überschallfeld. Hierbei ist die unmittelbare. 
Nachbarschaft der kritischen Geschwindigkeit der Behandlung mit der Charakeristikenmethode?) 
nicht mehr zugänglich. Besonders bei der Konstruktion rotationssymmetrischer Überschallstrahlen 
macht sich dies geltend. Deshalb wurde hierfür eine analytische Untersuchung angestellt. Dabei 
zeigte es sich, daß es außer den Strömungen, bei denen aus dem Parallelstrahl mit der kritischen 
Geschwindigkeit ein Überschallstrahl entsteht, auch solche gibt, die nach Unterschall führen. Bei 
= diesem zweiten Typ ist die Umkehrung, nämlich die Strömungen, die in Unterschall beginnend bei 
m. ° einem Parallelstrahl mit der kritischen Geschwindigkeit enden, besonders interessant. Dieser Pa- 
Br: rallelstrahl wird nämlich auf endlicher Länge gebildet, während exakte Unterschallparallelstrahlen 
auf endlicher Länge nicht zu realisieren sind. Um dieses Ergebnisses willen sei die einfache Rech- 
nung hier angegeben. Auf die Form des Ansatzes wird man bei der Betrachtung von Störungen in 
parallelen Überschallstrahlen geführt. Wir werden diese Untersuchung in unsere Rechnung ein- 
beziehen?), | 


u RN? m 


2. Grundgleichungen, Ansatz für die Behandlung von Störungen. 


Es seien x und y die Lagekoordinaten, bei rotationssymmetrischen Strömungen ist die x-Achse 
N die Achse der Rotationssymmetrie, 
& und 7 ein weiteres rechtwinkliges Koordinatensystem bei dem die 7-Achse mit einer 
Machschen Welle des Parallelstrahls, der Nullpunkt mit dem Nullpunkt des 
x, y-Systems zusammenfällt (vgl. Bild 1), 
der Geschwindigkeitsbetrag 
die örtliche Schallgeschwindigkeit, 
= arcsin a/w der Machsche Winkel, 
das Geschwindigkeitspotential, 
eine Zahl, die für ebene Strömungen null, fürrotationssymmetrische Strömungen 
eins ist, 
x das Verhältnis der spezifischen Wärmen bei konstantem Druck und konstantem 
Volumen, und 
I eine Vergleichslänge. 
Größen des Ausgangsparallelstrahls tragen einen Index 0. 
!) L. Prandtlund A. Busemann: Näherungsverfahren zur zeichnerischen Bestimmung von ebenen Strö- 


mungen mit Überschallgeschwindigkeiten. Stodola-Festschrift, Zürich 1929, und R. Sauer: Theoretische Einführung in 
die Gasdynamik S.108ff., S.137ff. Berlin, Springer-Verlag 1943. 


2) Der Ausbreitung von Störungen im stationären Überschallstrahl ist die nichtstationäre Ausbreitung von Zylinder- 


oder Kugelwellen in ruhendem Gas analog. Sie wurde von W. Hantzsche und H. Wendt untersucht. Im wesent- 
lichen folgen wir ihrem Vorgehen. 


O2 8 


VE Bl  Aasana dla Nana 9 Bel A 1 DA Dal a a Da nn > a dann Ef) 0 ae ni 


1 heißt ee ER 
„(a 0) + wi... 


Nierbi En Ditferentistiönen. BR Andides angegeben Die Schallgeschwindigkeit a erhält mar 
i aus. aD B ernou . lischen ee ie 


EN 


Bild1. zy- und ?„-Koordinatensystem. 


Der un bestimmt sich aus 3 

= +d=-d+ D ae De A MATT 
Im 5 n-System lautet die Potentialgleichung | 
P, sin &, — PD; cos a, 


De (a Ze) 28, 2 a ei u = £ ee 
nsin ag — &C0S a, 
Für den Ausgangsparallelstrahl ist 
355 a EN ET VE (3a) 
oder 
E D,— w,: (Esin a, + N cos a.) Ka ee Desan EIN. 
Macht man & und 7 mit / dimensionslos, d.h. setzt man 
& N 
| E — TE 3 N T ee TE ee (4), 
so ist 3” 
Di wer DES N COS ea ge ac) 


Die Linie &=0 ist Machsche Welle. Es sollen nun Strömungen untersucht werden, die 
längs dieser Ma ch schen Welle eine Unstetigkeit besitzen. Für & < 0 wird das Geschwindigkeits- 
potential durch (3c) gegeben, für & > O mag zu dem Ausdruck (3c) ein Zusatzglied hinzutreten. Die 
Gestalt des Zusatzgliedes ergibt sich daraus, daß bei der Fortpflanzung einer Unstetigkeit längs 
einer Mach schen Welle der Charakter der Unstetigkeit wahrscheinlich ungeändert bleibt und nur 
die Stärke eine Änderung erfährt. Wir setzen daher für& >0 


d=w,' 1 (Esinao + mcoso + f(n) € u er re $); 


_ hierin bestimmt der Exponent m die Art der Unstetigkeit, der Verlauf von fihre Stärke. Längs der 
Verzweigungslinie & = 0 ist das Potential und seine Ableitung stetig. Daraus folgt m > 1. Später 
wird sich sogar ergeben m > 2. Gleichung (5) stellt eine Welle dar, die sich für positive 7 in der 
Strömungsrichtung von der Achse entfernt, für negative 7 der Achse (bzw. bei ebenen Strömungen, 
der Symmetrielinie des Strahls) nähert. 


; 3) Hierbei sind höhere Potenzen in & bereits weggelassen. Allgemeiner mußte das Zusatzglied lauten E” [f, (n).+ 
han)e+th(n)®re..). 


TEE 
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Die Ableitungen von fnach n werden durch Striche gekennzeichnet. Es ist dann 


= wo (sin ap + mE" fm) 


Di = m (m 1)ErTE Sa) 


®, = ws (005 00 + E"-f’ (n)) ._ 


der — a3" 03-2 msin fg: €” Be ) 
+ Größen mit &” und höheren Potenzen von £. 
| Wir setzen dies ein und beschränken uns dabei in jedem Glied auf die niedrigsten Potenzen 
in &. Dabei ist zu beachten, daß cos a, für &, = 5 verschwindet, so daß dann Glieder, die eine höhere 
Potenz von £& als die Glieder mit cos a, besitzen, von Bedeutung sein können. Man erhält 
— f2: m? (m 1) (# +1) *° - sin. —f’ (cos ag + &*-f)-Am- E*T'- sin, | 


sIn2 
sin? &g 


(6). 


(f’ sin. &* —m - €". f- cos a) =0 | 


+f"'- €" (sin 2, — cos? &) + ß* 


3. Unstetigkeiten in Überschallstrahlen. 
Bei Überschallstrahlen ist cos, + 0. Indem man die niedrigsten Potenzen von & berück- 
sichtigt und durch —f’:2 m £* "sin a, cos a, dividiert, erhält man aus (6) 


m (m —1)- («+1 Br y 
IE ' I: Lug Be ee, TRETEN (7). 
COS &y Be 
Ist m < 2, so überwiegt für kleine & das erste Glied, man erhält 
: = 0 
d.h. Unstetigkeiten dieses Charakters sind in der Strömung nicht möglich. Für m > 2 erhält man 
of+p-T=o. 
N 
Hieraus für ebene Strömungen (ß = 0) 
f 77 C, 
und für rotationssymmetrische Strömungen (BP =1). 
1 
Von 


C, und C',, wie überhaupt alle in dieser Arbeit auftretenden O;, sind Integrationskonstante. 
1 


Demnach gehteine Unstetigkeit, diesich in Richtung auf die Symmetrieachse zu fortpflanzt, wie n 2 
nach unendlich. 
Für m = 2 ergibt sich die Bernoullische Differentialgleichung 


’ y p} % =[° 1 
a en 
Für ebene Strömung ergibt sich als Lösung 
ERSTE 
a BR On nie) 


Ist f positiv, so wird durch die Unstetigkeit eine Geschwindigkeitssteigerung und eine Ab- 
senkung des Druckes hervorgebracht, bei negativen f hat man es mit einer Verdichtung zu tun. In 
diesem letzten Falle ist ©, + n negativ und erreicht mit wachsendem ») den Wert 0. f wird dann 
negativ unendlich. Die Verdichtungswellen steilen sich zu einem Verdichtungsstoß auf. Es ist auf- 


Sa > 


« 


a cosa ID See 
2a@+D +0 In 


EAN ara sich auch bei rotationssymmetrischer Strömung eine Verdichtungswelle ii. 


wachsendem 7 zum Verdichtungsstoß auf. Für Wellen, die nach der Strahlmitte zum verhältsich 


1 
‚in unmittelbarer Nachbarschaft der ae die Singularität wie n De 


4. Unstetigkeiten im Strahl mit der kritischen erden 


Bei der kritischen Geschwindigkeit ist cos &, = 0, sinayu=1. Das £&, n-Koordinatensystem | 


fällt mit dem x, y-System zusammen. Bei ee Division durch &” erhält man aus (6) 


PER T Mar a 5 Bw 
aa n . 
/ F 


au Glieder niedrigsten Grades in & 


+3 mm): ee 


An Gegesatz zu der entsprechenden Gleichung (7) für Überschallstrahlen erhält man hier eine | 
gewöhnliche Differentialgleichung 2. Ordnung. Damit ist der Verlauf der Unstetigkeit nicht mehr 


allein von der Störung abhängig, die man auf der einen Seite des Strahles anbringt, sondern es wird 
noch eine andere Randbedingung, in unserem Falle die Störung auf der gegenüberliegenden Seite 
des Strahls, wirksam. Für den rotationssymmetrischen Strahl läßt sich die zweite Randbedingung 


am besten als Symmetrieforderung, f’= 0, für die Mitte des Strahls, erfüllen. Daß eine Diffe-, 


rentialgleichung zweiter Ordnung entsteht, ist im Einklang mit der Tatsache, daß für die kritische 
Geschwindigkeit die beiden Birhlungen, die die Charakteristiken in be ch Malen besitzen, 
zusammenfallen. 


Für m <3 erhält man aus der letzten Gleichung 


RR 
Störungen dieses Charakters sind also unmöglich. Für m > 3 ergibt sich die Gleichung 
rn 8 
mit den Lösungen _ - 
RER -f=-6+6'n 


für ebene Strömung und 

2 | J=9%+0s lan 

für rotationssymmetrische Strömung. Wegen der Symmetrie in der Mitte ist (, = 0, =0. 
- Für m=3 erhält man die = 


B: 1804 1)- ae ee RS (8) 


mit der Bedingung f’=0 für „= 0. Wir De zunächst U lolgendes, Ist F (£) eine Lösung der 
Differentialgleichung 


at 1dF 
EN ee Re (9) 
BrZE 
I dF 
mit dE —ı0 
für Z = 0, so sind die gewünschten Lösungen von (8) gegeben durch 
1 | 
— (2. F(C.: 
a 


Durch geeignete Wahl der früher eingeführten Vergleichslänge I läßt sich stets erreichen, daß 
die Integrationskonstante C,=1ist. Es gibt für F zwei verschiedene Gestalten F, und F,, die 


zusammen mit - in den Bildern 2 bis 5 dargestellt werden. Um einen Einblick in die Bedeutung 


dieser Lösungen zu erhalten, sei der ebene Fall genauer untersucht. 
13 


g, daß sich diese Brei nnie Dee für. m > 2 findet. Eine Unstetigkeit ie m—2 bedeutet 
_ einen Sprung in der Krümmung der Stromlinie. 


Bei a Ber Sapeisischer Strömung lautet die entsprechende Lösung Dr De 


> E R 
EN 


194 | Guderley, Störungen in Schall- und Überschallparallelstrahlen 


Dort erhält man als Lösung von (9) 


3 en FRE dT 4) 
=(5) :% Et 
z=+t1 


Unter Benutzung der oberen Vorzeichen läuft dann # zwischen einem positiven Wert, der sich 
aus T—= lerrechnet, und unendlich. £ variiert dabei zwischen 0 und einem endlichen Wert. In der 


0 02 04 06 08 10 0 02 04 06 08 70 
Bild 2. F, und 0,1 F, für ebene und rotations- Bild 3, Fi und 0,1 F, für ebene und rotations- 
symmetrische Strömung. symmetrische Strömung. 


A für rofationssym 


& für rolattonssymetrische 


elrische 
Strömung 


f für ebene Strömung 


0 22 04 06 08 10: =, 02 04 


06 98 10 


Bild4. F',für ebene und rotationssymmetrische Strömung. Bild. Fs für ebene und rotationssymmetrische Strömung 


“) Das Integral läßt sich auf ein elliptisches Integral zurückführen. Vgl. Jahnke-Emde ‚ Funktionentafeln. 


Leipzig 1933, S. 131. 
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hier Fortsetzungen über {= 1 hinaus denkbar. 


NEN PINAR 


mr 


Nähe der Stelle F= oo entspricht die Strömung dem Beginn der Me y erschen Strömung um die | 


EckeÖ), der zugehörige Wert £ entspricht der-Breite des engsten Querschnitts. Im übrigen ist F} 


stets positiv; diese Lösung stellt daher überall eine Überschallströmung dar. Mit den unteren Vor- 


zeichen variiert F, von einem negativen Wert über Null nach positiv unendlich: Interessant ist das 
Gebiet negativer Ro. Dieser Teil der Lösung stellt einen Unterschallstrahl dar, der in einen Strahl 
mit der kritischen Geschwindigkeit übergeht. Man muß dazu nur eine Stromlinie in geeigneter 


Weise als Rand ausführen. Benutzt man das gesamte negative Gebiet bis zur Stelle‘F', = 0, so ist 
; auf der Randstromlinie innerhalb unserer Näherung der Druck konstant. Diese Lösung läßt sich 


daher auch als Freistrahl realisieren. Damit ist gezeigt, daß ein Freistrahl mit der kritischen Ge- 


‚schwindigkeit sich auf endlicher Länge parallel richtet. Ganz entsprechend liegen die Verhältnisse 


bei Rotationssymmetrie. Die Integrationskonstanten wurden so gewählt, daß für den ersten Fall 
= 1bei F = oo im zweiten Fall bei F = Oliegt. Die Vergleichslänge I ist demnach einmal gleich 
der halben Breite bzw. dem halben Durchmesser des engsten Querschnitts, das zweite Mal gleich 
der halben Dicke des Freistrahls. N 3 


In den Bildern 2 und 3 ist für ebene und rotationssymmetrische Strömung F, und en so- 


dc 


wie, um einen Überblick über den gesamten Verlauf zu erhalten, 0,1-F, und 0,1- ee aufge- 


ß ar, : 2 2 
tragen. Im gewählten Maßstab ist für IE ein Unterschied zwischen ebener und rotationssym- 
metrischer Strömung nicht zu erkennen. Die Bilder 4 und 5 zeigen F', und ze Natürlich sind 


5. Zusammenfassung. 
= In Überschallstrahlen pflanzen sich Unstetigkeiten längs einer Mach schen Welle fort. Der 
Charakter der Unstetigkeit bleibt dabei erhalten, die Stärke ist aus einer Differentialgleichung erster 
Ordnung berechenbar. Dabei interessiert besonders das Verhalten der Unstetigkeit in der Nähe 
der Achse rotationssymmetrischer Strömung, außerdem läßt sich das Aufsteilen von Verdichtungs- 
wellen zu Verdichtungsstößen verfolgen. Auch in den Strahl mit der kritischen Geschwindigkeit 
hinein können sich Unstetigkeiten fortpflanzen. Die Stärke dieser Unstetigkeit wird hier durch 
eine Differentialgleichung zweiter Ordnung beherrscht. Von den beiden Lösungstypen stellt der 
eine die Expansion eines Schallstrahles nach Überschall hin, der zweite die Entwicklung eines Pa- _ 
rallelstrahls mit Schall aus einer Unterschallströmung dar. 


53) Th. Meyer: Diss. Göttingen 1908. 
Eingegangen im August 1944. 


Zusammenfassender Bericht. 


Über die Charakteristikenverfahren der Hydrodynamik. 


Teil 1... 
Von Kl. Oswatitsch in Graz. 
- (Aus dem Kaiser-Wilhelm-Institut für Strömungsforschung, Göttingen.) 
Das Gemeinsame der bisher bekannt gewordenen Oharakteristikenverfahren der Hydrodynamik wird 
= herausgearbeitet und ein besonders einfaches Verfahren angegeben für die praktischen Anwendungen. 
What they have in common the methods of characteristics of hydrodynamics being known hitherto is demon- 
sirated and an especially simple method is given for practical use. 
Ce quil y ade commun entre les methodes de caracteristiques de la hydrodynamique, qui nous sont connues 
jusqu’& ce jour est demontre, et une methode particulierement simple est donnee pour la pratique. 
Oömee CTABIIHX 10 CHX IIOP H3BeCTHEIMH MeTOXOB O0PaßoTku ['unponuHaMHEH TOoNBepTaeTtch 
o6paÖdorke Hu OCODeHHO YIPOIMeHHEIÄa CIOCOO 0O6pabor Eu Öyner ob’ABIIeH MIA IPAKTHYECKOTO 
HpuMeHeHHn. 
1. Einleitung. 

In den letzten Jahren sind verschiedene Charakteristikenverfahren auf dem Gebiete der 
Hydrodynamik entwickelt worden (z.B. von Döring, Guderley, Sauer, Schultz-Grunow, 
Tollmien und Schäfer). Öfters hat ein und dasselbe Problem eine Bearbeitung von seiten 
mehrerer Autoren erfahren. Im folgenden soll nach einer Darstellung der Grundlagen aller dieser 
Verfahren eine besonders einfache rechnerisch-zeichnerische Methode entwickelt werden, die sich 
in jedem Falle anwenden läßt. Diese wird am Charakteristikenverfahren für achsensymmetrische 

13* 
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196. 
stationäre adiabatische Überschallströmung und an der instationären Stromfadens 


_ verfahren und einer übersichtlichen Darstellung ihrer Grundlagen dienen. Es wird sich im folgn- 
‘ den daher nicht vermeiden lassen, altbekannte und schon öfters veröffentlichte Erkenntnisse 


 Oswatitsch, Über die Charakteristikenverfahren de 


YR n 


veränderlichem Querschnitt: näher erläutert. ira De ee 
Unsere Arbeit soll vor allem einer Vereinfachung in der Anwendung von Charakteristiken- 


Er 


noch einmal wiederzugeben. Auch ist das hier abgeleitete Verfahren naturgemäß den Methoden 
anderer Verfasser nahe verwandt. "Sofern uns deren Arbeiten zugänglich sind, werden wir im 


einzelnen auf solche Ähnlichkeiten hinweisen. i 


ee - 
RL nr Se 


2. Grundgleichungen. DETTEER 
Beschränken wir uns auf zwei unabhängige Veränderliche, also entweder zwei Ortskoordinaten 


-oder eine Ortskoordinate und die Zeit, so ergibt sich die Möglichkeit, folgende praktisch wichtigen 


Probleme zu behandeln: Die stationäre ebene und die stationäre achsensymmetrische Überschall- 


: strömung, die instationäre eindimensionale, achsensymmetrische (Zylinderwellen) und kugel- 


symmetrische (Kugelwellen) Unter- und Überschallströmung. Bei der instationären eindimen- 
sionalen Strömung lassen sich dabei zwei Sonderfälle unterscheiden. Die einfache ‚‚instationäre 


- Rohrströmung‘‘ (ebene Wellen) und die „instationäre Stromfadenströmung‘, bei welcher ver- 


änderlicher Querschnitt des Kanals zugelassen wird, die Geschwindigkeit und die Zustandsgrößen 
aber über den Querschnitt des Kanals als konstant angenommen werden. Natürlich kann die 
Ausbreitung ebener Wellen- als einfacher Spezialfall der ‚‚instationären Stromfadenströmung‘“ 
angesehen werden. E 
“Wir haben es bei den ebengenannten Strömungen mit folgenden Differentialgleichungs- 
systemen zu tun: Mr: 


‚Stationäre ebene und stationäre achsensymmetrische Strömung [5]): 


2 2 { 2 : 
(5), u+[1-3),=0; =). Gasdyn. Gl. (Kontinuität) 


> SEND UST E (= p..d 4 
u, + 9% = ER Er =y m Drehung. 

Die eingeklammerten rechten Seiten der Gleichungen gelten für die achsensymmetrische 
Strömung. Im Falle ebener Strömung sind im allgemeinen die üblichen Bezeichnungen gewählt. 
c ist die örtliche Schallgeschwindigkeit, o die auf die Masseneinheit bezogene Entropie. Diese 
hängt lediglich von der Stromfunktion y ab. Bei achsensymmetrischer Strömung ist x die Achsen- 
richtung, y die Richtung quer zur Achse. Entsprechendes gilt für die Geschwindigkeitskomponenten 
u und v. Die Stromfunktion bei der achsensymmetrischen Strömung ist definiert durch uoy=y,; 
voy= —y, (bei der ebenen Strömung wie üblich durch yy = we; ya = — vo). Die par- 
tiellen Ableitungen sind durch Anfügen der unabhängigen Veränderlichen als unterer Index 


- gekennzeichnet. 


‚Instationäre adiabatische Strömung mit einer Ortskoordinate [4]: 


c2 5 
MRS le -  Euler- Gleichung. 


1 1 | 
Wr + e 0, t+—0a=f(e). Kontinuitätsbedingung. 


In der Euler-Gleichung läßt sich ja unter Annahme adiabatischer Strömung das Druck- | 
glied durch ein Dichteglied (Dichte = o) ersetzen. zist bei Kugelwellen (Zylinderwellen) der Kugel- 
radius (der Zylinderradius). Die Funktion f(x) hat folgende Bedeutung: 


Für Kugelwellen Ste) an 
für Te lerellen @)= Ge R 
x 
für ebene Wellen (instationäre Rohrströmung) Tla)=D | 
und für die instationäre Stromfadenströmung =, : 


!) Die Zahlen in eckigen Klammern beziehen sich auf das Schrifttumverzeichnis am Schluß der Arbeit 
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R wobei unter F der Stromfadenquerschnitt zu verstehen ist, der zeitlich unveränderlich angenommen 
_ wird. Wir können unsere folgenden Betrachtungen auch auf zeitlich veränderliche Stromfaden- 


querschnitte ausdehnen und auch Strömungen behandeln, bei welchen sich die Entropie von 
Massenteilchen zu Massenteilchen ändert. Dies spielt bei der Berechnung von Strömungen mit 
Verdichtungsstößen eine-Rolle. Die rechte Seite unseres Gleichungssystems wird dann etwas 


BR komplizierter. Prinzipiellläßt sich aber das Problem genau so behandeln wie stationäre Strömungen 


mit Verdichtungsstößen. 


’ 


Allen genannten Problemen liegt also ein Gleichungssystem von folgender Gestalt zugrunde 3 


019,4 0:9, + %h, ta hy=4A(v, 9,9; h); a) 

6192 + 829, + b;h, + bıh,= B(x,y,9,h); ’ y 

Wir haben es mit einem System von Differentialgleichungen erster Ordnung zu tun, welche in 

den Ableitungen linear sind. Die Koeffizienten a,, b, hängen dabei lediglich von den abhängigen 

Veränderlichen g und A, nicht aber von den unabhängigen Veränderlichen x und yab. Die Glieder 

A und B können sowohl Funktionen von x und y als auch von g und A sein. Es handelt sich also 
um ein spezielles quasilineares System von Differentialgleichungen erster Ordnung. 

Bei Strömungen mit Entropieänderung steht auf der rechten Gleichungsseite allerdings 


noch eine Größe, welche von der Stromfunktion abhängt. Die Stromfunktion erhält man aber 


erst mittels eines Integrals über die abhängigen und unabhängigen Veränderlichen. Dies bedeutet 


wohl eine kleine Komplikation aber keine prinzipielle Erschwerung ‘der Aufgabe, wie wir noch - 


sehen werden. 


Wir wollen noch eine ziemlich verbreitete Bezeichnung einführen. Durch die Größen g. 


und k werden bei unserem Problem die Strömungszustände beschrieben. Dies geschieht entweder 


- mit Hilfe zweier Geschwindigkeitskomponenten bei stationärer. Strömung mit zwei örtlichen 


Koordinaten oder mit Hilfe der Geschwindigkeit und einer thermodynamischen Zustandsgröße, 
der Dichte go, dem Druck, der absoluten Temperatur oder auch mit Hilfe der Schallgeschwindig- 
keit, welche ja lediglich von der Temperatur abhängt und den Vorteil hat, ebenfalls die Dimension 
einer Geschwindigkeit zu haben. Wir wollen daher g und h ganz allgemein als ‚‚Zustandsgrößen“ 
bezeichnen und, wenn wir irgendwelche Funktionen in Abhängigkeit von g und h als Veränderliche 
betrachten, von der „Zustandsebene‘‘ sprechen. x, y seien hingegen die Ortskoordinaten, ganz 
gleich, ob darunter tatsächlich zwei örtliche oder eine örtliche und eine zeitliche Größe verstanden 
wird. Die x, y-Ebene nennen wir Ortsplan oder Strömungsebene. 


3. Eigenschaften hyperbolischer Differentialgleichungen. 


Wir denken uns nun folgende Aufgabe gestellt. Die Zustandsgrößen g und h des Gleichungs- 
systems (1) seien auf einer Kurve £& (x, y) = £, gegeben (£, sei eine vorgegebene Konstante). _Ge- 
fragt ist nach den Werten von g und A in der Umgebung dieser Kurve. Wir können die Aufgabe 
so lösen [1], daß wir uns krummlinige Koordinaten eingeführt denken, deren eine Schar 
durch £ (x, y) = konst. gegeben sei und deren andere Schar n (z, 4) = konst. lauten möge. Dabei 
wird natürlich vorausgesetzt, daß jeder Punkt x, y nur einem Punkt £, n zugeordnet sei und um- 
gekehrt, was sich mathematisch durch die Ungleichung 


RER 
Na> Ny 
ausdrückt. Da g und h auf der Kurve & (x, y) = &, bekannt sind, sind auch die Ableitungen g 


und Ah längs dieser Kurve gegeben. Führen wir im Gleichungssystem (1) also & und n anstatt x 
und y als unabhängige Veränderliche ein, indem wir 


9,—= gIe$, -F In Na USW. 


setzen, so erhalten wir in 
(@72, fe a,8,) get (a3 &, ie G,8,) hz = A—(aın, H Gay) In — (N, =F QNy) Rn > | 
(b1&, + b2&,)9: + (ba&, + da&y) he = B—(bına + bon) In — (ba + bar) An; | 
gerade zwei Gleichungen für die unbekannten Ableitungen g: und he. Wir können also g und h 
auf einer benachbarten Kurve £& (z, y) = &, + d£ berechnen: 
Kar) + gl mdd; 
h(&,+dS,n)=hl&,n) the, mdE. 
Auf &(x,y) = &--dE können wir nun wieder in Kurvenrichtung differenzieren und erhalten 


so in gleicher Weise die Werte von g und A auf einer Kurve & (2,9) = &+2d£. In Schwierig- 
keit kommen wir dabei nur an den Enden der Kurve, da sich dort die Ableitungen längs der Kurve 


praktisch nicht genau bilden lassen. 


(8) 


EU NEIN RENTEN 


N RETTEN 


 Oswatitse 


Mit Hilfe der Kramerschen Regel können wir die Werte von ge und h; direkt ange 
 Bs ist: RER Er R = Ai en 


A—(a,13 + an )9n—(a Ne + any) ns; Et as, 
Beben by) 9n — (da Ne +bsny Ihn; daöz + bdadu 


= a6, +Q2&,;5 0, tar. ee 
| retten A 
a,&,+0s£8,; A— (011g + @yy) In —(Qg Ns + Cara) in | G a, Er 


0,6, ta; tar 
b, En +b28,; b; fan b, Er 


9: und h: lassen sich also überall dort-berechnen, wo die Nennerdeterminante nicht verschwindet. 


Wir wollen uns im folgenden aber gerade für jene ausgezeichneten Sonderfälle interessieren, a 
denen die Nennerdeterminante gleich Null ist. Wir nehmen also die Kurve & (x, y) = &ı nicht als 


gegeben an und berechnen von ihr ausgehend die Werte der Zustandsgrößen g: und h; auf einer 


Nachbarkurve, sondern fragen uns, welche Bedingungen eine Kurve & (z, y) = konst. erfüllen 
muß, damit wir die Berechnung der Zustandsgrößen auf einer benachbarten Kurve auf die “5 
gegebene Art nicht durchführen können. Die Bedingung erhalten wir durch Nullsetzen der 
Nennerdeterminante von Gl. (4), was zu folgender Gleichung führt: R 


(a,b, —azb,)&: + (aba —asbı + 056, —agb;) E, Ey + (aba —asb)Ey— 0 5); 


d El N a 
Wir erhalten für die Richtung der Kurve £ (x, y).= konst. Tr — = eine quadratische Glei- 
chung, also in jedem Punkt gerade zwei Lösungen. Sind diese tatsächlich reell, ist also 
"D2— (a,b, —a,bı +36, —a,b5,)? —4 (a,b, —azb,)- (a5, — a,b,)>0 er, - (6), 


so sprechen wir von einem hyperbolischen Gleichungssystem. (Unter D wollen wir dabei im fol- 
genden, um Eindeutigkeit zu erzielen, stets den absoluten Betrag der Diskriminante verstehen.) 
Bei der ersten Gruppe von Beispielen, den stationären Strömungen, ist dies bei Überschallgeschwin- 


digkeiten der Fall. Bei den instationären Strömungen, der zweiten Gruppe, haben wir es stets 


mit einem hyperbolischen Differentialgleichungssystem zu tun. Ist der obere Ausdruck kleiner 
als Null, so werden die Lösungen für die Kurvenrichtung komplex. Wir sprechen dann von einem 
elliptischen Gleichungssystem. Beim Verschwinden der Diskriminante spricht man von einem 
parabolischen Gleichungssystem. Wir sehen, daß lediglich die Koeffizienten a,, b;, nicht aber 
die Funktionen A, B der Gl. (1), für den Typus des Gleichungssystems entscheidend sind. 

Im Fall eines hyperbolischen Differentialgleichungssystems, mit welchem wir uns in dieser 
Arbeit ausschließlich beschäftigen wollen, haben wir also in jedem Punkt zwei ausgezeichnete 
Richtungen, bei denen die Nennerdeterminante von (4) verschwindet. Das gibt zwei Scharen 
von Kurven, welche die Eigenschaft haben, daß bei ihnen GI. (5) in jedem Punkt erfüllt ist. Diese 
ausgezeichneten Kurven nennt man Charakeristiken. Wir wollen im folgenden die beiden 
Kurvenscharen mit 

&(x,y)=konst. und n(x, y)= konst. 


bezeichnen und nehmen also nachträglich an, daß wir in den Gl. (3) und (4) bereits die Charak- 
teristiken als neue Koordinaten eingeführt hätten. Zu berechnen sind die Charakteristiken durch 


Integration der Gl. (5). Da aber die Koeffizienten dieser Gleichung von den Zustandsgrößen g _ 


und A abhängen und diese gerade ,als Funktionen von & und % gesucht werden, können wir die 
Integration natürlich nicht durchführen. 
Betrachten wir irgendeinen Punkt des Ortsplanes, in dem die Zustandsgrößen g, h einen 


stetigen Verlauf nehmen, so haben die Ableitungen g;z, 9y Ra, hy irgendwelche endlichen Werte 


und folglich auch die Ableitungen in irgendeiner anderen Richtung, z.B. ge und he. Da aber 
in Gl. (4) die Nennerdeterminante verschwindet, ist dies nur möglich, wenn gleichzeitig auch die 
Zählerdeterminante dieser Brüche gleich Null ist. In der Zählerdeterminante können wir unter 
Verwendung der Lösungen von Gl. (5): 
E20, dh eds Baba 
coty =— = -7 et 2 2 Te a 

4 & 2(0,6,—a,b,) ee: 
U A 0,5, —a,b, + mb —ab,+D 
N 2(a,6,—a,b3)‘ 
einige Reduktionen vornehmen. y, und y, sind die Winkel, welche die beiden durch den betrach- 
teten Punkt gehenden Charakteristiken mit der positiven Richtung der z-Achse einschließen, 


Er Teer (7b), 


3 Er 
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RN. STE yı und tg Ya also die Richtungskonstanten der Charakteristiken in der Strömungsebene. Außer- e: 
un; benützen wir folgende Abkürzungen, deren geometrische Bedeutung wir noch kennen lernen | 
werden: REEL ENDE Ih SR 

| 2% ee D N 
a ET a EN 
ER 2(a,b6,—a,b,) a 
EEE A BR BT 
Rob des ae a „oe WB DB), 
} ; ‚2(a,b,—a3b1) | z Be 
Wir erhalten dann durch-Nullsetzen der Zählerdeterminante in Gl. (4) für die Änderung der Zu- Br, 
..standsgrößen auf der Charakteristik & = konst. folgende Beziehung: _ $ x, Ber 
g (ee) AbaBass nd Ba; | BE 
- G s)E ide abge, Ba ee dba — 0 EISEN 2), 
f wobei s die längs der Charakteristik gemessene Bogenlänge bedeutet. Für die Differentialquotienten ge 
F ae wir dabei eine Bezeichnungsweise, wie sie ähnlich in der Thermodynamik üblich ist. 
3 heißt die Änderung von % mit der Bogenlänge auf der Kurve &=konst. Entsprechend 2 
j erhalten wir für die Änderung von g und A auf der Charakteristik n = konst. die Formel: 
= («9) Ab Ba. Ab, Ba, 
= e 3 2 de Tan a NE ande os ERNE (Ib). 
Re Die Charakteristiken haben demnach nicht nur die Eigenschaft, daß die Ableitungen der 


Zustandsgrößen, welche von den Charakteristiken wegführen, nicht bestimmbar sind, die Charak- 
teristiken haben außerdem die Eigenschaft, daß die Ableitungen der Zustandsgrößen g,h längs 7 
der Charakteristiken ganz bestimmte Bedingungen, nämlich die Gl. (9a) und (9b) erfüllen. Diese a 
Bedingungsgleichungen für die Ableitungen längs der Charakteristiken nennt man Verträg- Be 
lichkeitsbedingungen. Die eben besprochene Eigenschaft ist sehr bemerkenswert. 
Während wir g und A auf irgendeiner Kurve beliebig‘ vorgeben und dann die Werte auf einer be- 
nachbarten Kurve daraus bestimmen können, können wir auf einer Charakteristik weder die E 
Werte von g und A beliebig vorgeben, sondern diese müssen die Verträglichkeitsbedingungen a 
erfüllen, noch können wir, nachdem die Werte vorgegeben sind, weiterrechnen. a 

Nehmen wir nun an, die Zustandsgrößen in zwei benachbarten Punkten P, und P,, welche 7 
nicht aufein und derselben Charakteristik liegen, seien durch g,, h, und g,, h, gegeben, so können a 
wir die Werte von g und h in einem Punkt P berechnen, der mit P, und P, gemeinsam auf je einer 
Charakteristik liegt (Bild 1). Wir nehmen dabei an, daß wir die Abstände so klein gewählt haben, 
daß wir die Differentialquotienten in den Punkten P, und P, durch Differenzenquotienten er- 
setzen können. Typisch für die Charakteristikenmethode, die wir 
jetzt im Prinzip zeigen wollen, ist es dabei, daß wir die Zustands- 
größen nicht in einem durch die Koordinaten x, y vorgegebenen 
Ort berechnen, sondern den Punkt, in dem wir die Zustandsgrö- 
ßen haben wollen, dadurch festlegen, daß wir verlangen, er solle 
mit den Ausgangspunkten P, und P, je eine Charakteristik 
gemeinsam haben. Wir kennen also auch den Ort des Punktes P 
zunächst nicht, sondern müssen bei einer Charakteristikenmethode 
stets sowohl diesen wie auch die Zustandsgrößen selbst im 
gefragten Punkt bestimmen. Den Ort des Punktes P können wir 
sofort so bekommen, daß wir. mit Hilfe von Gl. (7) die Richtun- 
gen der Charakteristiken in den Punkten P, und P, bestimmen 

- und die zwei Charakteristiken, von denen die eine zur Schar € = konst., die andere zur Schar 

n= konst. gehört, zum Schnitt bringen. Auf diese Weise erhalten wir von P, und P, ausgehend 
zwei Punkte. Beide können wir berechnen. Jenen, welcher mit P, gemeinsam auf der Charak- 
teristik & = konst. liegt, wollen wir mit P, bezeichnen (Bild 1). _ 

Für die Werte von g und h im Punkte P, erhalten wir folgende Gleichungen; wenn wir mit 
$s— 5; und s3s—s, die Längen der Strecken P,P, und ?, P, bezeichnen: 


Bild 1. Gitterpunkte 
in der Strömungsebene. 
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we hı+ os (s—s)=h,+ (2), —8;) 
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Wir haben hier also zwei Gleichungen für folgende Ableitungen im Punkte 2 
und die Ableitungen im Punkte, P;: wi ; we . Die Verträglichkeitsbedingung ( 


| 9a) ‚liefeı 


- eine weitere Bedingung für die beiden erstgenannten Ableitungen; entsprechendes gilt für (9b) e 
und die beiden letztgenannten Ableitungen, so daß sich im ganzen also vier Gleichungen mit 
vier Unbekannten ergeben und daraus die vier gefragten Differentialquotienten berechnet werden ; | 


können. Damit lassen sich aber auch die Zustandsgrößen im Punkte P, ermitteln. 


; "Sind also die Werte von g und % auf einer Kurve vor- K 
gegeben, so kann auf diese Weise das ganze Feld,welchsde 


von den Endpunkten A und B (Bild 2) ausgehenden Cha- 
rakteristiken einschließen, berechnet werden. Wir brauchen 
auf der Kurve nur jeweils von benachbarten Punkten auszu- 
gehen und in der früher beschriebenen Weise die Werte in 
neuen Punkten auszurechnen. Diese können wir dann wie- 
der zum Ausgangspunkt der Fortsetzung unserer Rechnung 
23: EFT machen. > | 
Ser ee Das eben beschriebene Berechnungsverfahren nennen wir 
system. ° in üblicher Weise ‚„‚Gitterpunktverfahren‘. Bei ihm werden die 
Werte der Zustandsgrößen in’ganz bestimmten Punkten der 


Strömungsebene als bekannt angenommen und daraus die Werte in neu zu bestimmenden Punkten . 
ausgerechnet. Wir haben es hier mit einer Verallgemeinerung der numerischen Integration einer | 


gewöhnlichen Differentialgleichung erster Ordnung zu tun. 

. „Dem Gitterpunktverfahren kann man das sogenannte „Felderverfahren‘‘ gegenüberstellen. 
Hier wird vorausgesetzt, daß die Zustände in ganzen Feldern konstant sind, wobei man sich die 
Felder durch Charakteristiken begrenzt denkt. Aus der Kenntnis der Größen g und hin bestimmten 
vorgegebenen Feldern kann man die Zustandsgrößen in neuen Feldern bestimmen. Als typisches 
Beispiel für diese Methode seidas Prandtl-Busemann- Verfahren erwähnt. Bei Berech- 


nungsmethoden, bei welchen aber Bogenlängen vorkommen, wie dies im folgenden im allgemeinen, 


der Fall ist, ziehen wir das Gitterpunktverfahren vor, da man beim Felderverfahren gezwungen 
ist, Mittelpunktsdefinitionen für die einzelnen Felder einzuführen, um den ‚„‚Abstand‘‘ eines Feldes 
vom anderen zu erhalten. Die von uns gezeigte Berechnungsmethode läßt sich im übrigen genau 
so für Felderverfahren durchführen. : | 
... Bei der Behandlung des achsensymmetrischen stationären Problems werden wir gleich- 

zeitig das ebene stationäre Problem als Gitterpunktverfahren behandeln und auf diese Weise 
das Gegenstück zum Prandtl-Busemann- Verfahren gewinnen. 

; Wir können uns nun fragen, unter welchen Umständen sich die Gl. (9a) und (b) ohne wei- 
= integrieren lassen. Abhängig von den Ortskoordinaten der Strömungsebene sind in diesen 
: leichungen lediglich die beiden Funktionen A, B und ferner die Bogenlänge auf den Charak- 
won s. Daraus ersehen wir, daß sich die genannten Gleichungen nicht etwa dann integrieren 
: a wenn A und B nur von den Zustandsgrößen g, R und nicht von x, y abhängig sind. Eine 

ntegration ist vielmehr nur dann möglich, wenn beide rechte Gleichungsseiten von (9a) und 
(9b) verschwinden. Dies liefert uns zwei Bedingungen für A und B, nämlich 


A(b,+b,cot y))+ B(a3+ a,coty,)—=0 \ 
A(b;+b,cot y,)+ B(a3 + a, cot y)—=0 ) 


Dieses Gleichungssystem hat nur dann eine nichttriviale Lösung, wenn die Koeffizientendeter- 
minante verschwindet, was nur eintritt, wenn 


(11). 


cot yy=coty, 


ist, d.h. wenn die beiden Charakteristikenricht i 
ungen zusammenfallen. Dies hab | i 
Gl. (2) ausgeschlossen. (11) hat also nur die triviale Lösung. DD 


4d=B=0 1. ms en ee 


Diese Bedingung ist notwendig und hinreichend für die Integri ] i 
grierbarkeit der Gl. (9) ohne Kenntni 
nn eleinpuseen und Anfangsbedingungen in der Strömungsebene. San, im Falle Ri 
en a ee stationären Überschallströmung und der instationären verdichtungs- 
ne nn 0) ae verwirklicht. Diese beiden Fälle lassen sich besonders einfach behandeln. 
ee = an einem in den Differentialquotienten homogenen linearen Differential- 
nt a ne homogenen Differentialgleichungssystem erster 
lediglich von den ae a an) 


a 
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Ds ee a rongpeyekurn Ben wir Ennehse, RE Die Methode i ist hier 
besonders einfach, und wir können sie - weitgehend für die Lösungsmethode beim inhomogenen 


er ausnützen. 


4. (liataktöristikenverfahren für das homogene quasilineare 
Differentialgleichungssystem. 


erschrndet die rechte Gleichungsseite von (9a) und (9b), so erhälten wir folgende zwei 
Gleichungen für die Änderungen von g und h längs der Charakteristiken, welche wir für den Spezial- 
fall des homogenen Gleichungssystems (1) (A= B=0) mit 


E—=konst. und n= konst.- 


bezeichnen wollen. Für diese gelten genau dieselben Formeln wie für £= konst. und n=konst, 


nur wird überall A= B=0 gesetzt. Die Gln. (9) bekommen dann die Gestalt: 


dh dh | N (138, DR 
() tea (= Ge Dt (d13a,b); 


tg 6, und tg ö, sind also die Richtungskonstanten der Charakteristiken des homogenen Systems $ 


in der g, A-Ebene. Sie hängen nach Gl. (8a, b) voraussetzungsgemäß lediglich von den Zustands- 
größen g und A ab. Wir können also (13a, b), abgesehen von mathematischen Schwierigkeiten, 


integrieren und ersehen daraus, daß wir beim homogenen Gleichungssystem (A= B=0) fol-. 


gende interessante Eigenschaft vorfinden.: Denken wir uns zwei Strömungsebenen. In beiden 
greifen wir Punkte heraus, bei welchen die gleichen Zustände g,, Ah, herrschen mögen und ziehen 
durch diese Punkte die Charakteristik & = konst. Wir suchen nun auf dieser Charakteristik in 
beiden Strömungsebenen Punkte auf, die im Wert von g übereinstimmen. Sie müssen dann auch 
den gleichen Wert Ah haben. 

Als einfaches Beispiel können wir bei der ebenen stationären adiabatischen Überschall- 


strömung den Überschallwirbel und die Überschallquelle betrachten. Als Zustandsgrößen wählen . 


wir Geschwindigkeitsbetrag und Geschwindigkeitsrichtung. Greifen wir nun jeweils den Punkt 


heraus; wo der Geschwindigkeitsbetrag gerade gleich der Schallgeschwindigkeit ist und der Winkel 


der Geschwindigkeit mitder x-Achse 9 = ist und suchen wir auf der durch diesen Punkt gehenden 
Charakteristik & = konst. jenen Punkt auf, in welchem der Geschwindigkeitsbetrag w gerade die 
Maximalgeschwindigkeit erreicht, so stimmen in diesem Punkt die Geschwindigkeitsrichtungen 
in beiden Strömungsebenen überein. 

Wir wollen uns nun die Frage stellen, wie ln das Bild der Charakteristiken £, 2) des 
homogenen Systems Gl. (7) in der Zustandsebene aussieht. Um dieses zu finden, müssen wir die 
Brüche £,/&z und n,/n« in &,/&, und »/ns überführen. Wie man durch Differenzieren leicht 
zeigen kann, ist 


& _ 
3 Iy—Iaz ie 
N te ee 
Z h,—h, in eb 
R © 


Setzen wir die Gl. (7) in (14) ein, so erhalten wir unter Berücksichtigung der Gl. (1) und (12) nach 
einigen Umformungen folgende Formeln für die Bilder der Charakteristiken des homogenen Systems 
in der g, -Zustandsebene: 


= 2(a4,b, —a,b,) RE . (15) 
_ m ub—mbtnb—agb+D_ 52) — cotö e 
ng 2 (0,6; —a,bı) dh/n 


Dies sind zwei gewöhnliche Differentialgleichungen für die Funktionen & (g, h) und 7 (g, h), die wir 
ohne Kenntnis der Anfangsbedingungen in der Strömungsebene in der Zustandsebene integrieren 
können, da auf der rechten Gleichungsseite lediglich Funktionen von g und %$ stehen. 

Wenn also das Differentialgleichungssystem (1) homogen ist, und die Koeffizienten a, b; 
wie vorausgesetzt lediglich von den Zustandsgrößen g, h abhängen, so liegen im Zustandsraum 
die Bilder der Charakteristiken der’Strömungsebene ein für allemal fest. Sie sind von den speziellen 
Anfangs- und Randbedingungen des Problems unabhängig. 

Wir wollen uns nun fragen, ob diese Bilder der Charakteristiken der Strömungsebene in 
der Zustandsebene auch dort wieder Charakteristiken sind. Zu diesem Zwecke leiten wir uns die 
Charakteristiken der Zustandsebene ab. 


4 


Hydrodynamik 


h Oswatit sc 


A hen SEEN FT ee ee 
OR * Anstatt mit den beiden Funktionen g (x, 9), h (&, y) in der Strömungse 3 : er 
Br wollen wir ‘jetzt das Funktionenpaar x (g, h) und y(g, h) in der Zustandsebene betrac .- Dat. 


können das homogene Gleichungssystem (1) für diese beiden Funktionen leicht transformieren. I 
‘ Durch Differentiation erhalten wir: Sr LERNEN 2 en 2 
u Mean Bar Me en 

= — (92 hy — Rz 9y) Yo: hy = (9: hy — Rz 9y) %u ER i 


Be. ed a e . ; 2 EEE ee 8 
0° Dies in Gl. (1) für A= B= 0 eingesetzt ergibt: BUT a € 
Re .- N ee 2 ES = ö (17); 

Ei A £ RT bb ytbin=l R Bi 

4 also ein homogenes lineares Differentialgleichungssystem erster Ordnung. Hätten wir ds 
Sam inhomogene Gleichungssystem (1) genommen, so hätten wir keineswegs ein inhomogen lineares 

‚a Gleichungssystem für & (g, h), y(g, h) erhalten, sondern die Funktionen A und B wären dann 

Su mit (9,%,—h,g,) multipliziert gewesen. Wir hätten also auch Glieder zweiten Grades in den 
Be Ableitungen bekommen. 


Wesentlich ist, daß wir in (17) ein richtig lineares Differentialgleichungssystem und kein i 
quasilineares vor uns haben. Denn die Koeffizienten ai, b; hängen nun ja voraussetzungsgemäß 
nur vön den jetzt unabhängigen Veränderlichen g und A ab. Eine Voraussetzung für die Trans- 


2 

AS EER formation ist übrigens auch die, daß man sich auf zwei unabhängige und zwei abhängige Ver- 
a % - änderliche beschränkt. Tun wir dies aber, so können wir die Eulersche Gleichung und die 
Be Kontinuitätsbedingung stets in diese Form bringen. BE ur 
e*: Die Charakteristiken des neuen Gleichungssystems (17) können wir nun sehr einfach gewinnen, 
E: wir brauchen lediglich in der Charakteristikengleichung (5) diejenigen Vertauschungen durch- 
rn 


zuführen, welche die homogenen Gl. (1) in die Gl. (17) überführen: Es ist also x durch 9%Y durch 
h, a, durch a, usw. zu ersetzen. Bezeichnen wir beide Charakteristikenscharen zunächst mit 
€ (g, h) = konst., so erhalten wir: Z 
0,5, + 35)&4+ (mdı— a6, + 06— a0), rt bit a b)ä—t ne (18) 

mit den Lösungen: 3 
ud dr een 
& 2 (a,5,— a,b,) 

Es ist leicht zu zeigen, daß die Diskriminanten von Gl. (5) und Gl. (18) übereinstimmen, wes- 
_ ‚halb wir in beiden Fällen denselben Buchstaben D verwenden können. i ’ 

Vergleichen wir die beiden Wurzeln von (18) mit den Gl. (15), so finden wir, daß die Glei- 
chungen paarweise miteinander übereinstimmen. Damit ist gezeigt, daß die Bilder der Charak- 
teristiken der Strömungsebene in der Zustandsebene auch in dieser wieder Charakteristiken sind. 
Es ist also & (x, y) = konst. dasselbe wie € (x, y) = konst. 

Weiß man das, so findet man die Bilder der Charakteristiken in der Zustandsebene auf 
dem eben begangenen Weg sehr leicht. Allerdings sieht man auf diese Weise nicht, wie die beiden 
Scharen in der Zustandsebene den beiden Scharen in der Strömungsebene zugeordnet sind. Da 
wir diese Zuordnung aber noch brauchen werden, haben wir beim Aufsuchen des Bildes der Charak- 
teristiken in der Zustandsebene nicht die Charakteristikengleichung (5) selbst in die Zustandsebene 
transformiert, was etwas einfacher gewesen wäre, sondern dies für jede einzelne Wurzel durch- 
geführt. : 

Nachdem wir einmal erkannt haben, daß die Charakteristiken bei unserer vereinfachten 
Aufgabe A= B=0 im Zustandsraum unabhängig von den Anfangsbedingungen festliegen, ist 
es leicht, einfache Verfahren zur Berechnung der Strömungsebene-anzugeben. Unsere Aufgabe ist 
also die, den Ort des Punktes P,, der mit P, und P, auf je einer Charakteristik gemeinsam liegt 
(Bild 1) und die Werte der Zustandsgrößen im Punkte P, zu bestimmen, wenn Ort und Zustand 
von P, und P, gegeben sind. 5 

Zunächst wollen wir uns mit einem Verfahren erster Ordnun g begnügen, 
indem wir die Richtung der Verbindungsgeraden P,P, mit der Richtung der Charakteristik im 
Punkte P, gleichsetzen. Entsprechendes machen wir auch für die Charakteristik Pe 

Die Richtung der Charakteristik im Punkte P, hängt nach’Formel (7) lediglich von den 
Zustandsgrößen g, k im Punkte P, ab, die ja bekannt sind. Wir brauchen also entweder eine 
Tabelle, aus der wir bei vorgegebenen g und h die Richtung der Charakteristiken ablesen können 
oder wir haben in der g- und A-Ebene ein Kurvensystem, daß die Eigenschaft hat, in jedem Punkte 
9, h die zugehörigen Richtungen nach Gl. (7a) und (7b) anzugeben. In beiden Fällen braucht die 
Arbeit, die entweder zur Ausrechnung der Tabelle oder zur Berechnung des Kurvensystems nötig 
ist, für jedes Differentialgleichungssystem nur ein einziges Maldurchgeführt zuwerden. Einedritte 


ar Möglichkeit, die Richtung einfach zu bestimmen, besteht darin, daß wir anstatt der Zustandsgrößen \ 7 . R 
‘ gund Airgendwelche Funktionen von g und h als neue Zustandsgrößen wählen, diesoeinfachsind, 


' können. Wählen wir z.B. als neue Zustandsgrößen 


einfach 


‚und » die Richtungskonstante einfach finden. 


RR Eh } N at BER DNz Bus a EP RE u 
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daß wir aus ihnen durch elementare Rechenoperationen die Richtungskonstanten sofort ableiten 
m=m(g,h) und n—n(g, RE, Nr ee Re A 
 hb—ab,+ ba ae D Ss N 


’ Mm — - und N NE R 2 
(0b, —aub,) 2(0b, —asb,). \ ae Bi. 


so können wir durch Bildung von Summe und Differenz der beiden neuen Zustandsgrößen m. ee 


Überhaupt besteht in der Tatsache, daß wir im allgemeinen gar nicht gebunden sind, g 


R und h als Zustandsgrößen beizubehalten,. eine Möglichkeit, die Richtungen der Charakteristiken _ 
"in der Strömungsebene einfach zu bestimmen. Wir können neue abhängige Veränderliche ein>r 


führen, die wir nicht von vornherein festlegen, sondern dadurch gewinnen, daß wir an sie gewisse 
Forderungen stellen. Wir können das Differentialgleichungssystem (1) mit A=B=0fürg 
und A ohne-weiteres in ein Differentialgleichungssystem für m und » verwandeln, dessen Eigen- 
schaften sich von jenen des Differentialgleichungssystems im wesentlichen nicht unterscheiden. 
Suchen wir die Charakteristiken in der m, n-Ebene, so finden wir für diese nach längerer Rechnung: 
2% _ agb, —aubz)n, + (0,5, —ab,+%b—a,b)m +Dm 
En  2lasba —asb)m, + (asd1 —a ds + %b5—a;b)m, + Dm, "|. (23) 
Nm > 2 (@3 b,—a,b)n,+ (a,b1—a, b,+ %b; —a;b,)n, —Dn, - 
Nm 2 (ad, —aub;)m, + (aıdı —a,bı+ a,b; —a;b,)m, —Dm, u 
Man kann nun z.B. fordern, daß die Charakteristiken in der Strömungsebene- stets senkrecht 
zur Charakteristik der anderen Schar im entsprechenden Punkt der Zustandsebene stehen 
(eine Forderung, die beim gewöhnlichen Prandtl-Busemann--Verfahren erfüllt ist), 
also: 


2 
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Setzen wir in (22) Gl. (21) und (7), letztere mit &=Eundn=pein, so erhalten wir nach einigen 
Reduktionen folgendes Differentialgleichungssystem für die neuen. abhängigen Veränderlichen 
m und nals Orthogonalitätsbedingung: 


(a,b, —asb,) m, — (ad; —a, b,)m, + (9b, —a,b,)n,—=0 (23) 


(a,b, — as bs) m, + (a,b, —az b,) m, + (@; bu — a,b), =0 KEN \ { E 


- Alle Lösungen dieses Differentialgleichungssystems geben uns im Zustandsraum Charakteristiken, 


welche die Eigenschaft haben, daß die Tangente im Punkte m, n an eine Charakteristik senkrecht 
steht auf der Tangente derjenigen Charakteristik des Strömungsraumes bei den Zustandswerten j 
m und n, welche der anderen Schar angehören. 

Wir sehen also, daß es eine ganze Anzahl solcher Lösungen gibt und wir können uns unter 
den Lösungen jene aussuchen, welche uns am praktischsten erscheint. 

Ebensogut wie die Forderung der Orthogonalität können wir eine andere Forderung stellen. . 
Am nächstliegenden ist es wohl zu verlangen, daß die Charakteristik in der Strömungsebene 
parallel zuihrem eigenen Bild im Zustandsraume verläuft. In diesem Falle ist: 


En N 
Das gibt nach einigen Reduktionen folgende Bedingungen für den Parallelismus: 
— (0,6, —a,b,) m, + (ad, — ad), + (ab, — a d)m—0; | (25) 
—(@b4 —a; 5) m, + (@1d5 — a, d1)N, + (Qı dz —a,b,)N,—0 ) ee : 
Alle anderen Bedingungen, die man sich noch aufgestellt denken kann, kommen wohl weniger 
in Frage. Esist für die praktische Rechnung eine naheliegende Forderung, daß die Charakteristiken 
in der’Strömungsebene ungefähr die gleiche Richtung haben wie ihre Bilder in der Zustandsebene. 
Dies erleichtert das Auffinden der zugehörigen Charakteristik. Daher fordert man auch, daß die 
Richtung der Charakteristiken parallel zum eigenen Bild des Zustandsraumes oder 
senkrecht zum Bild der anderen Schar im Zustandsraum verläuft. 
Wir haben verschiedene Möglichkeiten kennengelernt, den Ort des Punktes P, zu finden, 
wenn wir die Zustandsgrößen in den Punkten P, und P, kennen. Wollen wir nun die Zustands- 


e 
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Zustandsdiagramm, in welchem wir uns in gewissen Abständen die Charakteristiken ein für alle- 


_ mal eingezeichnet denken, aufzusuchen. Schreiten wir vom Bildpunkt P, längs der Charakteristik ei 
. F= konst. und vom Bildpunkt P, längs der Charakteristik 7 — konst. fort, so gelangen wir bei. BL 
richtiger Wahl der Fortschreitungsrichtung zum Bildpunkt von P,. Seine Koordinaten sind die 


gesuchten Zustandsgrößen. | 


Unter Umständen kann man durch richtige Wahl der abhängigen Veränderlichen das Ver- 


fahren so vereinfachen, daß auf das Zustandsdiagramm überhaupt verzichtet werden kann. Wählen 


wir z. B. bei der instationären adiabatischen Rohrströmung die Geschwindigkeit und die Schall- 


geschwindigkeit als abhängige Veränderliche, so werden die Charakteristiken im Zustandsraum 
Geraden. Wir brauchen hier im allgemeinen das Zustandsdiagramm gar nicht erst aufzusuchen, 
sondern können den Zustand des Punktes P, durch einfache Regeln aus den Zuständen von P, 
und P, berechnen, in dem wir im Kopf machen was sonst auf dem Diagramm ausgeführt wird. 


_ Wir wollen im folgenden aber in jedem Fallbei unserer geometrischen Darstellung mit Betrach- 
tung des Zustandsdiagramms bleiben und überlassen es dem Leser, gegebenenfalls von dieser 


geometrischen Behandlung abzurücken. 


Als recht praktisch erweist es sich vielfach, die Konstanten der Charakteristiken des homo- 


genen Systems E 
&(g, h)=konst.; n7(g, h)= konst. 


selbst als neue Veränderliche zu wählen. Geht durch P, die Charakteristik & (g, k) =£, und durch 
P, die Charakteristik 7 (9, k) =7j,, so wissen wir sofort, daß P, die Zustandswerte £, und 7, hat. 
Kennen wir diese, so können wir aus dem mit den Charakteristiken £ und 7) versehenen g, h-Zu- 
standsdiagramm sofort die gewünschten Zustandsgrößen entnehmen. Dies auf das Felderverfahren 
bei ebener Strömung angewendet, gibt uns gerade das Prandtl-Busemann- Verfahren. 

Nach der Behandlung des Verfahrens erster Ordnung für das homogene Gleichungssystem 
ist es nun leicht, ein Verfahren zweiter Ordnung anzugeben. Haben wir die Zu- 
standsgrößen im Punkte P, und P,, so erhalten wir die Zustandsgrößen im Punkte P, innerhalb 
der Zeichengenauigkeit des Zustandsdiagramms genau, auch wenn sich die Zustände beliebig 
stark voneinander unterscheiden. Das Verfahren der Auffindung der Zustandsgrößen g, und A, 


ist also auch beim Verfahren: erster Ordnung innerhalb der Zeichengenauigkeit des Zustands- - 


diagramms exakt. Verfeinerungsbedürftig ist aber die Bestimmung der Charakteristikenrichtung. 
Es ist einleuchtend, daß wir zu besseren Resultaten kommen müssen, wenn wir zur Festlegung 
des Punktes P, in der Strömungsebene die Charakteristikenrichtung nicht in den Punkten ?, 
und P, nehmen, sondern in Zwischenpunkten der Strecke P, P, und P, P,. Diese Richtung tragen 
wir dann von den Punkten P, und P, aus auf und nehmen im Schnittpunkt P, an.“ Im allgemeinen 
bildet ja die Charakteristik, auf welcher P, und P; liegt, einen flachen Bogen. Die Verbindung 
P,P; ist die Sehne dieses Bogens und wir setzen diese einfach richtungsgleich der Tangente des 
Bogens in einem Zwischenpunkt, wodurch im allgemeinen eine höhere Genauigkeit erzielt wird. 
Als Zwischenpunkte wählen wir die Halbierungspunkte der Strecken P,P, und P,P, im Zu- 
standsdiagramm. 

Beim Verfahren zweiter Ordnung suchen wir also zunächst die Punkte P, und P, im Zu- 
standsdiagramm auf. Der Schnittpunkt der Charakteristiken im Zustandsdiagramm gibt uns 
innerhalb der Zeichengenauigkeit exakt die Zustandsgrößen im Punkte P;. Um den Ort des 
Punktes P, zu finden, nehmen wir im Zustandsdiagramm die Halbierungspunkte der Strecken 
P, P,und P,P, und bilden in diesen beiden Punkten die Charakteristikenrichtungen in der Strö- 
mungsebene. Diese Richtungen von P, und P, aus aufgetragen, geben uns den Ort des Punktes P,. 

Wir haben nun also Verfahren erster und zweiter Ordnung zur zeichnerischen Lösung eines 
hyperbolischen Differentialgleichungssystems (1) kennengelernt. Vorausgesetzt war hierbei, daß 
es sich um Charakteristiken im Zustandsraum handelt, welche unabhängig von den Anfangs- 
und Randbedingungen der einzelnen Aufgabe sind, was gleichbedeutend mit der Forderung 
A=B=0 ist. Die hier geschilderten Methoden sind im wesentlichen schon längere Zeit bekannt 
und veröffentlicht. Wir sehen, daß es beim homogenen System kaum lohnt, ein Verfahren erster 
Ordnung anzuwenden, da das Verfahren zweiter Ordnung nicht nennenswert komplizierter ist. 


5. Charakteristikenverfahren für das inhomogene quasilineare Differential- 
gleichungssystem. 


Den Ort des Punktes P, bei einem Verfahren erster Ordnung beim inhomogenen 
Gleichungssystem (1) (wenigstens eine von den beiden Funktionen A, B darf nicht überall ver- 
schwinden) zu finden, ist nicht schwieriger als früher. Denn die Richtung der Charakteristiken 
hängt hier nach Gl. (7) genau so wie beim homogenen Gleichungssystem lediglich von den Zu- 
ständen ab. Wir können die Charakteristikenrichtungen in den Punkten P, und P, also genau so 
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. liegt. &2,85,und „9, sind dabei beim Verfahren 
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bestimmen wie früher, etwa dadurch, daß wir auf Grund der Bedingungen der Orthogonalität 


oder des Parallelismus die Charakteristikenrichtungen einfach aus dem Zustandsdiagramm, in 


welchem wir uns die Charakteristiken des homogenen Gleichungssystems eingetragen denken, 
entnehmen. | 


Anders ist es aber mit der Bestimmung der Zustandsgrößen im Punkte P,. g und h müssen 
sich längs der Charakteristiken nach Gl. (9a) und (9b) ändern. Werden diese gewöhnlichen Dif- 
 ferentialgleichungen als Differenzengleichungen geschrieben, so erhalten wir, wenn wir uns P, 


wieder auf &= konst..und P, auf. n = konst. gemeinsam mit P, liegend denken: 


Ab,—Ba .  4Ab,—Ba = 
Er ee BD my _ Au Ba | BEN e 
2 hı—tg ö, (gg — 91) = a,b, ın Yı Abs —arb; c08 Yı (85 s)=S, (26a), 
n R ER Ab,—Ba, . Ab,— Ba, N | Er or n 
art h,—tg 0%, (g; — 9) = Be dub, Tr 2b, —a,b, cos 9] (3; —%)= 8, (26 b). 


. Die Größen ö, und ö, hängen dabei lediglich von den Zustandsgrößen g und h, die Ausdrücke 


S, und 8, außerdem noch von den Ortskoordinaten ab. Bei einem Verfahren erster Ordnung 
müssen wir uns dabei tgö, und den Klammerausdruck in (26a) im Punkte. P, gebildet denken. 
Entsprechendes gilt für Gl. (26b) und den Punkt P,;. 

. Um nun die Zustandsgrößen im Punkte P, zu finden, wollen wir den Gl. (26a) und (26b) 
durch Einführen neuer Koordinaten eine einfachere Form geben. Wir machen dies am besten 
so, daß wir die neuen Koordinaten in den gewöhnlichen Differentialgleichungen (9a) und (9b) 
einführen und erst dann wieder zu Differenzengleichungen übergehen. Zunächst wollen wir als 
unabhängige Veränderliche die Charakteristiken des homogenen Gleichungssystems & und 7 
wählen und erhalten nach längerer Rechnung dann folgende Differenzengleichungen 


für ge konsiet euer IT re 


h fürn konto mn ee en el bi 
Für das homogene Gleichungssystem 8, = 8, = 0 g 

führt dies zum trivialen Resultat, daß &, und &, 
gleich sein müssen, d h. ?P, und P, liegen ge- 
meinsam auf der Charakteristik & (x, y) = &. Ent- 
sprechendes gilt für P, und P,. Handelt es sich 
aber um ein inhomogenes Gleichungssystem, so 
sagen die beiden Gleichungen aus, daß der Punkt 
P, auf den Kurven 


E@h)=&4+&8, und nl h)=n+mS: 


erster Ordnung im Punkte P, und P, zu nehmen, 
es sind also Konstante. Wir finden also P, genau 
so wie beim homogenen Gleichungssystem, nur 
dürfen wir nicht von den Punkten P, (£,, 7) und 
P; (£,7s) längs der Charakteristiken = konst. 
und n7= konst. fortschreiten, sondern von zwei __ : 

korrigierten Punkten P/(&+ 58,7) und """ "nhomogenen Ditferentialgleichungssystem. 
Pi(&m+mS,), Bild3. © > 

Die Ausdrücke 


=, 25 [4, = Bo, Au—Ba,. | 

&,8, = An: 08 ER Re (s:—5}) ; 
Ab,— Ba, Ab,— Ba, . | 

1.8, = cos EIHASL (33 — 8) 


hängen bei einem Verfahren erster Ordnung von folgenden Größen ab: 1. von den Zustandsgrößen 
in den Punkten P, und P,. Mit ihrer Hilfe können wir &,, 27, @, 5; und y; berechnen. 2. Von 
den Ortskoordinaten der Punkte P, und P,. Nimmt man diese noch dazu, so kann man auch 
A und B ermitteln. 3. Von den Bogenlängen der Charakteristikenstücke P,P,;und P,P;, s33;—s, und 
83—8;. Wählen wir die Maschenweite so klein, daß wir die Bogenelemente durch Geradenstücke 
ersetzen können, so können wir s,—s, und s,—s, dem Ortsplan ohne weiteres entnehmen, wenn 
wir. vorerst den Ort des Punktes P, bestimmt haben. Wählen wir die Charakteristiken des homo- 


genen Systems selbst als Zustandsgrößen, so können wir die Gl. (27a) und (27b) als Gleichungen - 


zur Ermittlung des Zustandes in P, auffassen und diesen direkt in das Zustandsdiagramm ein- 
tragen. Dieses Vorgehen findet man vor allem bei Guderley. 
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Br Der Vorgang bei einem Verfahren erster Ordnung zur Lösung des inhomogenen Gleichungs- ” 
Be. systems ist also folgender. Zuerst bestimmen wir den Ort des Punktes P; genau auf dieselbe E 
RER Weise wie beim Verfahren erster Ordnung für das homogene Gleichungssystem unter Verwendung “re E 
0 der Charakteristiken des homogenen Gleichungssystems. Dann bestimmen wir die Abstände = 
Br. $;—s| und s,—s, im Ortsplan und können nun mit Hilfe dieser Abstände, der Ortskoordinaten 


- und der Zustandsgrößen der Punkte P, und P, die Korrekturen &,8, und nn S, berechnen. Jetzt 
suchen wir die korrigierten Punkte P/ (&, + &,8,,71,) und P, (&,72+ 7 8,) im Zustandsdiagramm _ 
auf und bestimmen von diesen Punkten ausgehend, im übrigen aber genau so wie beim Verfahren 
für das homogene Gleichungssystem, den Punkt P,. Wir haben hiermit ein Lösungsverfahren 
erster Ordnung für alle inhomogenen Differentialgleichungssysteme von der Form der GEL: 

er Wie wir sehen werden, ist die hier allgemein geschilderte Methode im Einzelfall recht einfach zu 

N: handhaben. i EUER 

Br Für die Praxis kommt es vor allem darauf an, die Korrekturen möglichst einfach zu be- 

Br: rechnen, was im wesentlichen auf eine einfache Berechenbarkeit der Klammerausdrücke (26a) 

und (26b) herausläuft. Wir spalten diese daher möglichst in Teile auf, welche lediglich von den 

B: Zustandsgrößen g und %, und in Teile, welche lediglich von den Ortskoordinaten abhängen, auf. 

h Die Abhängigkeit von den Ortskoordinaten ist bei den von uns am Anfang genannten Beispielen 

Be in der Regel recht einfach. Etwas komplizierter ist hingegen die Abhängigkeit von den Zustands- 

fL.- größen. Hier können wir uns aber dadurch helfen, daß wir uns die erforderlichen Funktionen ein 

& für allemal berechnen! Wie dabei vorzugehen ist, entscheidet man am besten beim einzelnen 

Verfahren. Ari 
Ebensogut wie in der g, A-Ebene können wir die Rechnungen natürlich auch in der m, n- 
Zustandsebene durchführen. Dabei denken wir uns m und n so gewählt, daß die Orthogonalitäts- 
bedingung oder die Bedingung des Parallelismus erfüllt ist. Die Berechnung der Korrektur &,8, 
und 7,8, erfolgt in derselben Weise wie früher; denn wenn die Zustandsgrößen m und n bekannt 
sind, sind ja auch die Zustandsgrößen g und h gegeben. Die Charakteristiken des homogenen 
> Systems & und n.sind in der m, n-Ebene dann natürlich gegenüber ihrem Verlauf in der g, k-Ebene 

. ....° in irgendeiner Weise verzerrt. | 

Wir haben zur Ermittlung des Punktes P, in der Zustandsebene den Punkt P, längs einer 
Charakteristik der einen Schar nach P/ und den Punkt P,längs einer Charakteristik der anderen 
Schar nach P/ verrückt. Es ist aber im allgemeinen gar nicht notwendig und auch nicht immer 
am einfachsten, die Verschiebungen der Punkte längs Charakteristiken durchzuführen. Uns 
stehen hier noch andere Möglichkeiten offen, wie wir im folgenden zeigen wollen. 

Denken wir uns, tgö hänge nicht von A sondern lediglich von g allein ab. Diese Annahme 
ist nicht sehr speziell, denn wir können uns die abhängigen Veränderlichen g und h so gewählt 
denken, daß sie die gewünschte Eigenschaft erfüllen. Außerdem kann diese Wahl der Veränder- 
lichen erfolgen ohne die Orthogonalitätsbedingung oder eine andere Bedingung zu verletzen. 
Denn wir wollen g und A gar nicht als Koordinaten im Zustandsraum verwenden, als solche wählen . 
wir etwa m und n, sondern wollen im Zustandsraum lediglich die Linien g = konst. und h = konst. 

. herausgegriffen denken. Diese bilden in der m, n-Ebene zwei Scharen von irgendwelchen krumm- 
linigen Kurven. Hängt tgö lediglich von g ab, so können wir die Gl. (26a, b) so deuten, daß wir 
uns denken, es wäre beim homogenen Gleichungssystem der Punkt P/ mit den Koordinaten 
9, und A, + $, und der Punkt P} mit den Koordinaten g, und A,-+ S, gegeben und der Punkt 
P; ist zu suchen. Wir finden ihn wieder, indem wir vom Punkt P/ -auf der Kurve E= konst. 
und vom Punkt P, auf der Kurve 7= konst. fortschreiten. Im, Schnittpunkt finden wir dann 
den Punkt P.,. 


Bild5. Strömungsebene mit Stromlinien. 


Bild 4 (links). Eine zweite Möglichkeit zur Bestimmung 
des Zustandes im Punkte Pi: 


Die Methode bleibt genau dieselbe wie früher, nur werden die Punkte P. | 
» > und P, nun 
anderer Richtung verschoben, um zu den Punkten P; und P; zu gelangen (Bild 4). Diesen Beten 
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letzten Punkten kommt weiter keine geometrische Bedeutung zu, sie sind lediglich Rechengrößen. 
Man sieht dies allein schon daraus, daß diese Verschiebungen lediglich den Zweck haben, Korrek- 
turpunkte P; und P, zu finden, welche auf den durch den Punkt P, gehenden Kurven &=&, 
und n=n,liegen. An welchen Stellen dieser Kurven die Punkte P/ und P/ dabei liegen, ist 


lediglich aus Erwägungen der Zweckmäßigkeit gegeben. Die Korrekturen sollen eben möglichst 
einfach berechenbar sein. 


Bei stationären und instationären Strömungen mit Entropieänderungen können Größen. 
vorkommen, welche in komplizierterer Weise von den Ortskoordinaten abhängen. Bei der statio- 
nären Strömung mit Verdichtungdstoß z. B. kommt die Änderung der Entropie mit der Strom- _ 
funktion vor, eine Größe, welche von der Stromfunktion,selbst und mithin vom Stromlinienverlauf 
abhängt. Wir sind also genötigt, Stromlinien bei der Konstruktion mitzuzeichnen. Denken wir 
uns die Stromlinien y=y; und y= y; +, etwa bis zur Verbindungsstrecke von P, und P, durch- 
konstruiert (Bild 5). Dies geht, weil ja die Werte bis in diesen Teil der Strömungsebene als bekannt 


vorausgesetzt sind. Wir können nun durch Interpolation > ‚ welches ja.nur von w abhängt, 
in P, und P, bestimmen und haben dann alle zur Beechime Kastieke Werte. Beim Verfahren 
do 
dy 

Wir wollen uns nun noch ein Verfahren zweiter Ordnung zur Bestimmung 
von Ort und Zustand des Punktes P, bei inhomogenem Gleichungssystem verschaffen. Hier sind 
verschiedene Wege denkbar. Wir können z.B. zunächst Ort und Zustand von P, mit Hilfe des 
Verfahrens erster Ordnung bestimmen. Nennen wir den so gewonnenen Punkt P/. Nun bestimmen 
wir den Ort des Punktes P, genauer, indem wir die Charakteristikenrichtungen in den Halbierungs- 
punkten der Strecken P,P/ und P, P; im Zustandsdiagramm nehmen. Dies gibt uns innerhalb 
des Verfahrens zweiter Ordnung bereits den richtigen Punkt P,. Jetzt bilden wir mit Hilfe des 
Punktes P, noch einmal s,— s, und s,— s,, ermitteln die örtlichen Koordinaten an den Halbierungs- 
punkten von P,P, und P, P,in der Störungsebene und die Zustandsgrößen in den Halbierungs- 
punkten von P, P/ und P, P/ und berechnen nun die Korrekturen 8, und S, noch einmal. Mit 
Hilfe dieser genauen Korrekturen finden wir nun die endgültigen Zustandswerte von P;. 


Im allgemeinen zeigt es sich, daß die Korrekturen in großen Gebieten unserer Rechnungen 
ziemlich unbedeutend sind. D.h. die Verrückungen der Punkte P, nach P/ und P, nach P/ 
im Zustandsdiagramm sind gegenüber dem Abstand der Punkte P, und P, im.Zustandsdiagramm 
nur klein. In diesem Falle ist eine andere Iteration vorzuziehen, welche schneller zum Ziele führt. 
Wir nehmen als ersten Schritt nicht das Verfahren erster Ordnung des inhomogenen Gleichungs- 
systems, sondern das Verfahren zweiter Ordnung des homogenen Gleichungssystems. 

Wir tragen die Punkte P, und P, im 
Zustandsplan ein und suchen den Punkt auf, 
den wir beim homogenen Gleichungssystem 
erhalten würden. Diesen nennen wir nun 
P, und nehmen ihn als erste Näherung der 
Lösung des inhomogenen Gleichungssystems. 
Die weiteren Schritte sind nun ähnlich wie 
früher. In den Halbierungspunkten der 
Strecken P,P; und P,P; des Zustandsdia- 
Sramms ermitteln wir die Charakteristiken- Bild 6. Verfahren zweiter Ordnung zur Bestimmung des Ortes 
richtungen, tragen diese im Ortsplan von P} und Zustandes von P, beim homogenen Differentialgleichungs- 
und P,ausgehend auf und finden so den Ort nz 
des Punktes P/ (Bild 6). Bis hierher ist alles 
eine Wiederholung des Verfahrens zweiter Ordnung für das homogene System. In den Halbierungs- 
punkten der Strecken P, Pi und P, Pi des Ortsplanes und der Zustandsebene bilden wir nun die 
Korrekturen £,$, und 7,8, und kommen durch Verschieben längs der Charakteristiken zu den 
Korrekturpunkten P/ und P/ in der Zustandsebene. Den Punkt P, im Zustandsdiagramm 
finden wir nun wieder ausgehend von den Punkten P/ und P/ durch Fortschreiten längs der 
Charakteristiken des homogenen Systems. Liegen die Punkte P, und P, im Zustandsdiagramm 
nicht unmittelbar nebeneinander, so müssen wir nun noch die Charakteristikenrichtungen in den 
Halbierungspunkten der Strecken P,P, und P,P, der Zustandsebene bilden und finden so den 
Punkt P, im Ortsplan. Eine abermalige Korrektur ist hier überflüssig, da wir bei dieser Methode 
die. Verschiebungen als verhältnismäßig klein zu den Abständen der Punkte P, und P, im Zu- 
standsdiagramm vorausgesetzt haben. Ist diese Voraussetzung nicht erfüllt, so verwenden. wir 
besser das erstgeschilderte Verfahren zweiter Ordnung. 

Natürlich könnten wir als erste Näherung für das Verfahren zweiter Ordnung des inhomo- 


erster Ordnung ist es ja nur folgerichtig, auch in den Punkten P, und P, selbst zu nehmen. 


- aber deshalb nicht, weil das Verfahren zweiter Ordnung 
Rechenaufwand erfordert, aber den Vorteil hat, uns ö 
wie beim Verfahren erster Ordnung kann auch beim Verfahren zweiter Ordnung die Korrektur 
in der Zustandsebene durch Verschieben der Punkte P, und P, längs anderer Kurven als den 
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genen Systems auch das Verfahren erster Ordnung des homogenen Systems nehmen. Dies lohnt 


Charakteristiken (Bild 4) erfolgen. Dies werden wir uns in den Anwendungen zunutze machen. 
Auch das Vorkommen von Größen, welche von der Stromfunktion abhängen, macht keine 


| - zusätzlichen Schwierigkeiten. Wir extrapolieren die Stromlinien aus dem bekannten Gebiet 


heraus einfach ein Stück weiter (Bild 5), so daß wir die Stromfunktion an den Zwischenpunkten 


von PÄ—P; und P,—P, bestimmen können. 


Eingegangen am 2. 7.1944. 


NACHRICHTEN 


Th. Pöschl 65 Jahre alt. 

Dr. Theodor Pöschl, ordentlicher Professor für 
Mechanik und angewandte Mathematik an der Tech- 
nischen Hochschule Karlsruhe, wurde am 6. September 
dieses Jahres 65 Jahre alt. 

Gebürtig aus Graz, studierte Pöschl vornehmlich an 
der Hochschule dieser Stadt unter dem maßgeblichen 
Einfluß von F. Wittenbauer, dem Schöpfer der ‚‚Gra- 
phischen Dynamik“. Diesem Pro lemkreis gehörenauch 


einige der frühen Arbeiten Pöschlsan. Die Verbunden-. 


heit mit seinem Lehrer zeigte sich zudem darin, daß 
Pöschl später die Wittenbauersche Sammlung von, ,‚Auf- 
gaben aus der Technischen Mechanik“ neu herausgab. 

Seit 1916 war Pöschl Ordinarius an der Deutschen 
Technischen Hochschule in-Prag. Aus der Lehrtätig- 
keit an dieser Hochschule erwuchs sein ‚Lehrbuch der 
Technischen wfechanik‘‘, das ihn bald weithin bekannt 
machte. Im Jahre 1928 folgte Pöschl einem Ruf an 
die Technische Hochschule Karlsruhe. Hier setzte er 
seine fruchtbare -Lehrtätigkeit fort und entwickelte 
eine bedeutende literarische Produktion. Neben einer 


‘großen Anzahl von Zeitschriftenaufsätzen aus allen 


Gebieten der Mechanik entstanden die Artikelim Hand- 
buch der Physik über den Stoß und über die techni- 
schen Anwendungen der Stereomechanik, eine Neu- 
auflage der ‚‚Technischen Mechanik“, dazu deren 
2. Band, der die Festigkeitslehre behandelt, ferner die 
schon erwähnte Bearbeitung der ‚‚Aufgaben aus der 
Technischen Mechanik‘ , die ‚‚Getriebelehre‘‘ und der 
Beitrag ‚‚Mechanik der nichtstarren Körper‘‘ in Müller- 
Pouillets Lehrbuch der Physik, der auch ins Englische 
übersetzt wurde. Der daran anschließende Bericht über 
„Mechanik der festen Körper‘ für „Die Physik in regel- 
mäßigen Berichten‘‘ markiert den Übergang zu einem 
neuen Arbeitsgebiet, der Physik der Werkstoffe; ihm 
widmete Pöschl im wesentlichen die Arbeit der Jahre 
nach 1937, in denen er von den, Nationalsozialisten 
seines Lehramts enthoben- war. 

Bei der Wiedereröffnung der Hochschule nach der 
Kapitulation übernahm Pöschl zunächst das Amt des 
Prorektors, im folgenden Jahre berief ihn das Ver- 
trauen seiner Kollegen zum Rektor. Mit unerschütter- 
ter und unerschütterlicher Arbeitskraft meistert der 
Jubilar die vielfachen Pflichten, die ihm als Lehrer einer 
riesigen Hörerschar, als Repräsentant der Hochschule 
und in manch anderem ehrenvollen aber auch arbeits- 
reichen Amt obliegen. Groß ist auch die Zahl der 
wissenschaftlichen Pläne und literarischen Entwürfe, 
die ihn bewegen. Wir dürfen noch manche ‚Frucht 
seiner ausgedehnten Kenntnisse, seiner großen Erfah- 
rung und seines nie erlahmenden Fleißes erwarten. 

Die unwandelbare Liebenswürdigkeit, Heiterkeit und 
Hilfsbereitschaft Pöschls wissen alle zu schätzen, die je- 
mals mitihm in persönliche Berührung gekommen sind. 


Karlsruhe, 7. September 1947. K. Klotter. 


Nach Berechnung des Punktes P, müssen die 
Stromlinien dann noch einmal genauer um das neu berechnete Stück weitergezeichnet werden. 
Der zweite Teil folgt in einem der nächsten Hefte. 


N 


Tagung in St. Louis, E 

In der Zeit vom 31. Juli bis 2. August d. J. fand in 
St. Louis eine Tagung über nichtstationäre Probleme 
der 'Gasdynamik statt. Außer deutschen Forschern 
waren Gelehrte aus den USA, aus Belgien, England, 
Frankreich und der Schweiz anwesend. Es wurden fol- 


"gende Vorträge gehalten: 


SCHARDIN (Weil/Rh.) — Einführung. 


SAUER (Weil/Rh.) — Mathematische Methoden zur 
Behandlung nichtstationärer Probleme der Gasdy- 


namik. i 
SAUTER (Weil/Rh.) — Streng lösbare eindimensionale 
Probleme. 

RICHTER (Weil/Rh.) — Verdünnungs- und Verdich- 
'tungsfächer. en 

HEINZ (Weil/Rh.) — Reflexion ebener und kugelsym- 
metrischer Druckwellen. 

SCHULTZ-GRUNOW (Aachen) — Theoretische und 
experimentelle Untersuchungen zur nichtstationären 
Rohrströmung. i - 

WEISE (Kiel) — Ausbreitung von Kügelwellen. 

WECKEN (Weil/Rh.) — Pseudostationäre Probleme. 

SAUTER (Weil/Rh.) — Über die Entropie und die 

“ Stabilität von Stoßwellen. 

WECKEN (Weil/Rh.) — Theorie der Mach-Welle. 

SAUER (Weil/Rh.) — Ausbreitung von Schallwellen 
in stationären Strömungsfeldern. 

BURKHARDT (Kiel) — Dissoziation und Ionisation 
von Gasen durch Steßwellen. 

SAUTER (Weil/Rh.) — Über das Leuchten von Stoß- 
wellen. 
KOHLER (Horb/N.) — Stoßwellen in Gasen mit an- 

geregten inneren Freiheitsgraden. 

SCHARDIN (Weil/Rh.) — Wirkung von Stoßwellen 
auf mechanische Objekte, 


WILLMS (Weil/Rh.) — Vermessung von Stoßwellen. - 


OERTEL (Weil/Rh.) — Dichtemessungen in Stoß- 
wellen mittels Ionensonde. _ 


BIERMANN (Hamburg) — Stoßwellen auf der Sonne. 


„Prof. Dr. L. FÖPPL wurde zum Rektor der Tech- 
nischen Hochschule München gewählt. 


Dr. H. SCHUBERT wurde zum a. o. Prof. der an- 
gewandten Mathematik an der Universität Rostock 
ernannt. 


Dem Leiter des Forschungsinstitutes für Mathematik 
und Mechanik in Leningrad Prof. SMIRNOW wurde 
wegen seiner Verdienste um die mathematische Wissen- 
schaft der Lenin-Orden verliehen. : > 


des homogenen Systems kaum mehr 
fters eine Iteration zu ersparen. Genau so 


Verantwortlich für den Inhalt: Prof. Dr. Fr. A. Willers, Dresden; für den Verla ı H.K S i 

Verlag: Akademie-Verlag GmbH., Berlin N 4, Chausseestr. 106, Fernsprecher: 42 50 01 Verlag ak ee 

App. 275), Postscheckkonto: Berlin 35 021. Bestell- und Verlagsnummer dieses Heftes: 1009/25— 277). Die Zeitschrift 

für angewandte Mathematik und Mechanik erscheint monatlich. Bezugspreis: vierteljährlich RM 15.— zuzüglich Bestell- 

geld, Einzelheft RM 6.—. Verantwortlich für den Anzeigenteil: Curt F. W. Schreiber, Berlin. Druck: 

Julius Beltz, Langensalza. (Th. B.-Nr. 04011/16b) — Nr. 1095. — Veröffentlicht unter der Lizenz-Nr. 245 der 
Sowjetischen Militär-Verwaltung in Deutschland. / 


